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Abstract 

This  report  Is  a  translation  from  the  Russian  of  a  survey 
article  by  L.D.  Faddeyev^  which  appeared  in  Uspekhi  Matem.  Nauk., 
]>,  No.  k   (88),  p.  57  (1959)-   Our  own  interest  in  this  article 
lies  in  its  relevance  to  the  inverse  scattering  problem,  -  that 
is,  the  problem  of  determining  information  about  a  medium  from 
which  an  electromagnetic  wave  is  reflected,  given  a  knowledge  of 
the  reflection  coefficient.   Similar  questions  concerning  scat- 
tering phenomena  in  other  branches  of  physics,  e.g.,  in  quantum 
mechanics,  can  be  investigated  by  means  of  the  same  theory.   We 
have  therefore  thought  it  worthwhile  to  reproduce  and  distribute 
the  translation.  A  good  indication  of  the  contents  is  given  in 
the  Introduction. 


Introduction 

This  paper  is  devoted  to  a  survey  of  the  following  fundamental  prob- 
lem arising  in  the  quantum  theory  of  scattering:  the  solution  of  the  equa- 
tion satisfying 

d^  2 

LiJ/  = g  ^ix,-k)   +   q(x)t(x,k)  =.k  v(x,k),  ■  (0.1) 

dx 

ijf(0,k)  =  0  (0.2) 

behaves  asymptotically  like 

\|/(x,k)  R=i  C(k)sin(kx  -  7i(k)),  (0.5) 

provided  the  potential  q(x)  decreases  sufficiently  fast  as  x  tends  to  in- 
finity; to  what  extent  does  assigning  the  function  Ti(k)  detemiine  q(x)  and 
how  are  they  related.   This  problem  is  one  example  of  the  general  question 
concerning  the  relationship  between  the  S-matrix  and  the  energy  operator 
in  scattering  theory.   The  operator  L  defined  by  equation  (O.l)  and  condi- 
tion (0.2)  is  the  simplest  example  of  the  energy  operator  occurring  in  scat- 
tering theory,  and  the  function  S(k)  =  e   ^^    ^   the  simplest  example  of 
the  S-matrix  or  scattering  operator. 

Introduced  for  the  first  time  by  Wheeler,  the  S-matrlx  has  since  been 
frequently  applied  in  scattering  theory  following  the  publication  of  Helsen- 
berg's  papers'- J.   in  these  articles,  the  following  time-independent  defini- 
tion of  the  S-matrlx  was  given.  A  wave  function,  describing  the  steady  state 
of  a  system  (for  simplicity,  we  restrict  ourselves  to  a  system  of  two  parti- 
cles), has  an  asymptotic  representation  In  the  space  variables 
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Y 


V-Ko 


^^  +  ^2  {O.k) 


where  r  is  the  distance  between  particles.  ^     and  Y   are  respectively 

outgoing  and  incoming  waves.   Therefore,,  1     contains  the  factor  e    and 

-ikr 
^2  ^^  proportional  to  e    ^  k  being  the  wave  number  characterizing  the 

energy  of  the  state  f .   The  quantity  relating  the  amplitudes  of  these  two 

functions  is  called  the  S-matrlx.   In  our  illustration^ 


^(x,k);^^[e^^-i^(^)-  e-i^+i^(^)' 


(0.5) 


l.e.^  the  first  term  corresponds  to  an  outgoing  wave^  the  second  to  an  in- 
coming wave^  and  their  amplitudes  aj:e  related  by  the  factor 

S(k)  =  e-2i^(^).  (0.6) 

This  expression  constitutes   the  S-matrlx  for  our  example. 

Heisenberg's   theory  of  the  S-matrlx  was   further  developed  in  the  papers 
of  MiZ^llerL   ^   who  gave   a  time-dependent  definition  of  the  S-matrlx  which 
physically  was  more  meaningful.      Since  then  this  time-dependent  formulation 
of  the   scattering  problem  has   received  a  great   deal  of  attention   (see   the 
appendix)    and  can  be  stated  in  the   following  way.      The   energy  operator  of 
a  system  consists   of  two  terms 

L  =  L     +  V,  (0.7) 

o  \        >  ' 

where  L  corresponds  to  the  energy  of  the  free  particles  and  V  to  the  inter- 
action energy.  Long  before  collision,  i.e.,  for  negatively  infinite  time,  the 
state  of  the  non-interacting  particles  is  described  by  a  vector  Q  (t)  whose 
dependence  on  time  is  determined  by  the  operator  L  : 


In  conformity  with  the  established  convention,  a  state  vector  will  be  some 
element  of  Hilbert  space  in  which  all  operators  act. 
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fl_(t)  =  e  ''"°  n_,  (0.8) 

n  being  a  constant  vector  characterizing  the  initial  state  of  the  system. 
For  finite  time,  the  state  n(t)  is  a  solution  of  the  Schroedinger  equation 

i  ^^  =  Ln(t)  =  (L^  +  V)n(t),  (0.9) 

and  is  req\alred  to  take  on  the  initial  state  Q.   (t)  in  the  sense: 

lim     lln(t)  -  n_(t)||  =  0.  (0.10) 

t  -►  -00 

Over  a  large  interval  of  time  after  collision,  the  particle  motion  again 
becomes  free  so  that  asymptotically  as  t  -►  <» 

||n(t)  -  n_^(t)||  -  0,  (0.11) 

-iL  t 

vhere  0  (t)  =  e   °  Q   .      The  manner  in  which  the  asymptotic  state  vector 
+  + 

changes  determines  the  nature  of  the  scattering  process.   The  operator  S 
which  relates  the  asymptotic  vectors  Q,     and  Q,       according  to  the  formula 

a     =  SCI  ,  •  (0.12) 

+     - 

is   called  the   scattering  operator  or  S-matrlx. 

In  Sec,    5,    it  will  be  shown  that  this   formulation  holds   for  the  example 
in  question,    the  S-matrix  being  given  by  the   function  S(k)   =   e        *^    ' 
occurring  in  the  time-independent  definition.      This   fact  typifies   a  general 
aspect  of  the  stationary  and  nonstatlonaiy  formulations   of  the  S-matrix  in 
case  both  definitions   are  valid. 
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Heisenberg  came  to   consider  the  S-matrix  as   a  means  of  overcoming  the 
difficulties   encountered  In      modem  relatlvistic  theory  of  elementary 

particles.     He  felt  it  vas  necessary  to  Introduce   a  new  fundamental   con- 
stant having  dimension  length.      Therefore^    he   analyzed  the   current  theory 
and  rejected  as   unobservahle   those   notions  which  contradicted  the   Idea  of 
a  fundamental  length.      Only  those   quantities   experimentally  observable  would 
be  put   in  a  future   theory.      In  this   sense^    the   S-matrix  satisfies   the   require- 
ments  of  Helsenberg.      It   depicts   the  wave   function   at   large   distances    and 
is  thus  not   contrary  to  the   fundamental  length  hypothesis.     Moreover,    the 
scattering   cross-sections,    which  can  be   directly  measured  experimentally, 
are   expressed  in  terms   of  the   elements   of  the   S-matrix.      Helsenberg  also 
conjectured  that  the  discrete  energy  levels   corresponding  to  bound  states 
of  the  particles   should  be  determined  by  the   analytic   continuation  of  the 
S-matrix  into  the   complex  energy  plane. 

Connected  with  Helsenberg 's   supposition  that  the  S-matrix  was  more 
fundamental  than  the  Hamlltonlan  was  the   question  of  clarifying  the   rela- 
tionship between  these  two   characterizations   of  a  system.      In  particular, 
in  what  sense   should  one  define  the  Hamlltonlan  on  the  basis   of  the  S-matrix 
when  both  of  these  notions  were  used  in  a  theory.      In  addition  to  the 
theoretical  aspect,    the  inverse  problem,    i.e.   the  reconstruction  of  the 
energy  operator  from  its   S-matrix,    could  be   of  great  practical  value  in 
the  interpretation  of  experimental  scattering  data  as  well  as   in  the 
determination  of  various  properties   of  the  particles  which  are  not  directly 
measurable. 

The  simplest  example  in  scattering  theory  is  the  radial  equation  which 
describes  the   scattering  by  a  fixed  spherically  symmetric   center: 


i_  ^(x,k)  ^i./'iL^Jj.   +  q(^))  ^(x,k)  =  k2^(x,k).    (0.15) 
dx  \   X 


This  reduces  to  the  case  already  mentioned  for  ^  =  0.   The  first  attempts  at 
solving  the  Inverse  problem  were  undertaken  by  Fr^dberg'-  -■  and  Hylleraas  i-  -^  . 
They  worked  out  a  formal  procedure  using  series  whose  convergence  was  very 
likely.   However_,  Bargmann'-  -■  constructed  explicit  examples  In  which  differ- 
ent potentials  gave  rise  to  the  same  S(k)  and  to  the  same  discrete  energy 
levels.   This  showed  that  a  potential  which  led  to  given  energy  levels  and 
scattering  functions  S(k)  could  not  be  reconstructed  In  a  unique  manner. 
LevlnsonL  -I  showed  that  this  lack  of  uniqueness  must  be  related  to  the  exist- 
ence of  a  discrete  spectrum.   To  wit,  he  proved  that  the  reconstructed  poten- 
tial Is  unique  when  there  Is  no  discrete  spectrum.   The  precise  mathematical 
reason  for  this  was  explained  by  Marchenko i-  -■  who  showed  that  the  S-functlon 
determines  the  continuous  portion  of  the  so-called  spectral  function  of  equa- 
tion (O.l).   To  find  the  spectral  function  when  there  Is  a  discrete  spectrum. 
It  Is  necessary  to  give  not  only  the  location  of  the  eigenvalues  but  also, 
for  example,  the  derivatives  of  the  corresponding  normalized  elgenfunctlons 
at  X  =  0.  Marchenko '- -'  also  showed  that  the  spectral  function  determines  a 
potential  uniquely.   Thus,  Marchenko  related  the  problem  In  question  to  the 
so-called  Inverse  Sturm-Llouvllle  problem  which  already  had  been  discussed 
In  the  mathematical  literature.   An  analogous  result  was  obtained  at  approxi- 
mately the  same  time  by  Borg^  J.   By  developing  Levlnson's  method,  Jost  and 
KohnL-'^-I^L  J  Independently  came  to  the  same  conclusion  concerning  the  reason 
for  this  lack  of  uniqueness.   In  addition  to  the  original  potential  q(x), 
they  gave  an  explicit  formula  for  a  family  of  potentials  having  the  same 

We  will  call  the  scattering  operator  in  our  example  the  S-function. 


-  6  - 


S-function  and  the  same  discrete  energy  levels.      An  analogous   formula 
was   obtained  by  Holmberg'-     -' . 

A  procedure  for  explicitly  constructing  a  potential  without  the 
slngulaj:-lty  l{l  +  l)/x     from  its  spectral  function  was   formulated  by 

ri2i 

Gelfand  and  Levitan'-     -^  .      They  reduced  the  problem  to  a  Unear  integral 
equation  and  gave  sufficient   conditions  in  terms  of  the   spectral  function 
guaranteeing  that  it  be  the   spectral  function  of  some  equation  with  a 
potential  from  a  given   class.      The   results   of  Gelfand   and  Levitan  on  the 
Inverse  Sturm-Llouville  problem  were  immediately  applied  to  the  inverse 
scattering  problem  by  Jost   and  KohnL     -^    and  Levlnson"-     -' .      In    [ij]  ^    a 
formula  was  given  for  a  family  of  equivalent  potentials  having  one   and 
the  same  S-functlon  and  discrete  energy  levels.      Stricter  conditions   (both 
necessary  and  sufficient)    on  the   spectral  function  were   obtained  by  Kreln"-     -■ 
His  paper  completed  the  general  problem  of  reconstructing  equation  (O.l) 
from  its   spectral  function.      However,    since  the  passage  from  the  S-function 
to  the  spectral  function  was  not  entirely  trivial,    there   still  remained 
the  unanswered  question  concerning  the  nature  of  the  set  of  possible  S-func- 
tions   corresponding  to  the  potentials   from  a  given  class.      This  problem 
was  solved  by  KreinL     ^    and  Marchenko'-     -'   who  gave   a  series   of  conditions 
in  terms  of  the  Fourier  transform  of  the   function  S(k)    -   1.     Marchenko 
showed  that  the  potential  q(x)   possesses   the  same  properties   for  x  tending 
to   zero   and  infinity  as  does  the  derivative  of  this  Fourier  transform. 
The   final  inequalities   deduced  by  Marchenko  permitted  him  to  formulate 
necessary  and  sufficient   conditions  on  the  S-function  ensuring,  that   a  poten- 
tial from  a  given  class  would  correspond  to  it. 
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After  the  basic  papers  of  Gelfand  and  Levitan^  Krein  and  Marchenko^ 
a  great  deal  of  work  vas  devoted  to  carrying  over  their  results  to  an 
equation  containing  the  singular  term  l{l  +  l)/x  ,    an  equation  over  the 
interval  -  »  <  x  <  «  ,  a  system  of  equations^  and  the  relativlstic  equa- 
tions.  A  brief  survey  is  given  in  the  appendix. 

It  IS  Interesting  to  note  that  in  the  USSR  the  inverse  problem  has  been 
studied  on  the  vhole  by  mathematicians  whereas  abroad  almost  exclusively  by 
physicists  who  merely  use  the  method  of  Gelfajid-Levitan  as  interpreted  by 
Levlnson,  Jost  and  Kohn.  An  explanation  of  the  general  features  of  this  me- 
thod, and  its  application  to  the  solution  of  various  problems,  was  undertaken 
in  a  series  of  papers  by  Kay  and  Moses  L  J  ~  L  J .   These  authors  used  the  gene- 
ral concept  X)f   transformation  operators  developed  by  FriedrichsL  -^ '  ^    -^ . 

Recently  work  has  been  devoted  to  applying  the  results  of  the  inverse 
problem  in  the  interpretation  of  experimental  scattering  data!-  -I  ■"  L  ^J  •>  L  J. 

Thus,  the  inverse  scattering  problem  for  the  simplest  example  of  the 
radial  equation  was  solved  during  the  course  of  about  a  decade^  and  a 
large  amount  of  literature  has  been  devoted  to  it.   In  the  present  survey, 
we  will  endeavor  to  give  the  results  of  a  majority  of  these  papers  in 
their  most  general  form.   This  will  make  clearer  the  way  in  which  the 
basic  results  are  carried  over  to  other  problems.   In  this,  the  general 
approach  to  transformation  operators  developed  by  Friedrichs  and  applied  to 
the  inverse  problem  by  Kay  and  Moses  will  play  an  essential  role. 

All  of  the  basic  results  in  the  inverse  scattering  problem  could  be 
obtained  using  one  of  the  methods  of  Gelf and-Levltan,  Marchenko  or  Krein. 
Our  presentation  will  not  stick  to  any  particular  one  of  these  methods,  but 


rather  at  different  points  vlll  make  use  of  different  methods.  We  will 
attempt  to  establish  their  connection  considering  that  each  one  of  these 
methods  explains  different  aspects  of  the  mathematical  structure  of  the 
entire  problem. 

Due  to  the  large  material  content^  all  mathematical  proofs  will  not  be 
carirLed  out  in  a  completely  rigorous  fashion.  Many  of  our  considerations 
will  be  of  a  heuristic  nature  where  the  justification  of  details  would  re- 
quire greater  backgroiind  than  in  other  more  standard  proofs.   We  will  never- 
theless use  these  heuristic  proofs  in  order  not  to  obscure  the  conceptual 
side  of  this  work  with  lengthy  mathematical  arguments.  We  are  confident 
that  the  physicist-reader  will  find  our  discussions  completely  convincing 
and  that  the  mathematician  will  be  able  to  reconstruct  the  deficient  proofs 
so  as  to  make  them  completely  rigorous.   On  the  other  hand^  we  have  tried 
to  state  theorems  in  their  most  precise  fonn. 

Let  us  briefly  outline  the  basic  ideas  and  plan  of  the  survey.   The 
first  13  sections   are  devoted  to  the  solution  of  the  inverse  scattering 

problem  for  the  operator  L  defined  by  equation  (O.IJ)  and.  the  condition 

o 
\|/(0)  =  0  without  the  singularity  l{l  +   l)/x  .   We  treat  the  operator  L  as 

a  perturbation  of  the  operator  L  defined  by  the  differential  expression 

2       2 
L  t  =  -  <3.  \(/(x)/dx  and  the  same  condition  t(o)  =  0.   According  to  Friedrlchs^ 

a  transformation  operator  U  is  defined  as  the  solution  of  the  operator  equa- 
tion 

LU  =  UL  .  (0.1^) 

o 

This   Implies   that  any  transformation  operator  which  has   an  inverse 
generates   a  similarity  transformation  of  the  perturbed  operator  into  the 
■"onpertiirbed: 
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U""'xU  =  L  .  (0.15) 

o 

This  transformation  operator  U  replaces  the  elgenf unctions  of  the  continu- 
oxis  spectrum  of  the  operator  L  in  all  considerations.  Roughly  speaking,  its 
kernel  is  obtained  by  expanding  the  eigenfunctions  of  the  continuous  spectrum 
of  L  in  terms  of  the  eigenfunctions  of  the  operator  L  . 

In  Sec.  h   and  Sec.  5^  it  Is  shown  that  such  transformation  operators 
exist  for  our  example  and  that  the  completeness  theorem  for  the  eigenfunc- 
tions of  L,  proven  in  Sec.  2,    can  be  written  down  in  terms  of  the  transforma- 
tion operator  in  the  form 

U¥U*  =  I,  (0.16) 

(for  simplicity,  we  have  restricted  ourselves  here  to  the  case  when  L  has  a 
discrete  spectrum;  in  the  text,  this  restriction  is  not  imposed) .   Here,  W 
is  a  positive-definite  self-adjoint  operator  commuting  with  L  .  W  determines 
the  'normalization'  of  the  corresponding  operator  U. 

A  characteristic  feature  of  our  example  is  that  eimong  the  transformation 
operators  there  exist  Volterra  operators  of  the  form 

X 

U^f(x)  =  (I  +  K)f(x)  =  f(x)  +   rK(x,y)f(y)dy.  (O.l?) 

o 

The  operator  W  corresponding  to  U„  is  constructed  using 

si 

^^^)  =  M(k)M(-k)  '  (0-18) 

where  M(k)  is  a  certain  function  introduced  In  Sec.  1.   One  might  call  M( v/x) 
the  determinant  of  the  operator  L  -  XL.  In  fact,  in  Sec.  2  it  is  shown  that 
this  function  appears  in  the  denominator  of  the  resolvent  kernel 

The  subscript  B  used  throughout  and  which  transliterates  into  English  as  V 
stands,  of  course,  for  Volterra. 
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of  L  axid  determines  its  singularities.  These  consist  of  a  branch  cut  and 
poles  corresponding  respectively  to  the  continuous  spectrum  and  points  of 
the  discrete  spectrum. 

In  Sec.  3j  it  is  shown  that  the  time -dependent  form.ulation  of  the 
scattering  problem  is  valid  for  our  example  provided  that  L  is  taken  as 
the  energy  operator  of  the  free  particles.   The  corresponding  scattering 
operator  is  determined  by  using  the  function 

SM  -  101  .  (0.19) 

In  Sec.   6,    it  is   shown  how  to  establish  the  relationship  between  W(k) 
and  S(k)   with  the  help  of  (0.l8)    and   (0.19). 

In  Sec.    8^    on  the  basis   of   (0.l6)    and  the  triangularity  of  the  kernel 
K(x,y),    a  linear  integral  equation  is  obtained  connecting  the  kernels   of 
W  and  K.      In  Sec.    9,    this   equation  is   studied   and  the   inverse  problem  is 
solved  for  the   case  when  L  doesn't  have   a  discrete   spectrum.      Supplementary 
facts  necessary  for  a  treatment  of  the  general  case   are   cited  in  Sec.   12. 

The   approach  described   corresponds   to  the  Gelfand-Levltan  method. 
Another  procedure,    related  to  Marchenko's  method,    is  based  on  the   applica- 
tion of  the  operator  V     =  I  +  A  introduced  in  Sees.   ^I-  and  7: 

B 

00 

Vgf(x)  =  f(x)  +  r  A(x,y)f(y)dy.  (0.20) 

X 

This  operator  is  not  a  transformation  operator  in  the  sense  of  the  general 
definition.  However,  its  relationship  to  the  transformation  operator 

U  =  U  W  =  (U*)"  is  established  in  Sec.  7.   From  this  relation  and  (0.I6), 

B    B      B 

it  follows  that  the  next  identity  holds  for  the  operator  V^: 
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Vg(l  -  F)V|  =  I,  (0.21) 

vhere  the  operator  F  can  be  constructed  directly  using  the  function  S(k). 
By  means  of  (0.2l),  a  linear  integral  equation  is  deduced  vhich  relates 
the  kernel  A(x,y)  to  the  function  S(k),  thus  permitting  one  to  solve  the 
inverse  problem.   This  integral  equation  is  used   in  Sec.  10  to  investigate 
the  connection  between  S(k)  and  q(x).   Several  aspects  of  Krein's  method 
are  illustrated  in  Sec.  11.   In  Sec.  15^  the  construction  of  some  opera- 
tor L  from  a  knovn  operator  L  when  the  S-fun'ction  of  L  differs  from  that 
of  L,  by  a  rational  factor  is  considered.   This  is  an  important  problem 
in  applications.   In  Sec.  ik   and  15^  the  results  deduced  are  carried  over 
to  the  radial  equation  (0.15)  when  t   >  0. 

In  order  not  to  Interrupt  the  presentation^  we  will  not  refer  to 
original  papers  in  the  text.   The  literature  is  cited  in  a  special  appen- 
dix.  A  series  of  comments  are  made  there  and  a  brief  sur^^ey  is  given  of 
work  done  on  the  inverse  scattering  problem  which  is  not  included  in  the 
text. 

For  her  valuable  hints  and  corrections,  the  author  expresses  his 
deepest  thanks  to  0.  A.  Ladi zhenskaya  who  read  through  the  survey  in  manu- 
script. He  would  also  like  to  thank  Z.  S.  Agranovich  and  V.  A.  Marchenko 
for  their  courtesy  in  allowing  him  to  look  at  their  monograph  prior  to  its 
publication. 
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!•    The  solutions  cp(x,s),  f(x^s)  and  their  relationship,  existence  and 
inequalities.   The  function  M(s)  and  its  properties. 

In  this  section,  some  basic  features  of  solutions  of  the  equation 

-y"  +  q.(x)y  =  s  y,   s  =  a  +  It,  (l.l) 

are  assembled  which  will  be  utilized  in  the  subsequent  presentation.  In 
all  lemmas,  it  is  assumed  without  further  mention  that  q(x)  is  a  locally 
summable  function  and  satisfies  the   condition 

00 

/     x|q(x)  |da  =   C  <  oo    .  (1.2) 

o 

The  solutions  q)(x,s)  and  f(x,  s)  are  determined  by  the  conditions: 

cp(x,s):   (p(0,s)  =  0,   q5'(0,s)  =  1,  (l.j) 

f(x,s):     lim  e"  ^^f(x,s)  =  1.  ,^   ^. 

Equation  (l.l)  and  the  conditions  (1.5)  and  (1.4-)  are  equivalent  to  the 
following  integral  equations: 


cp(x,s)=^HL^+   r   Bin  s(x--t)  ^(,)^(,^3),,^  (1.5) 


f(x,s)  =  e^^^  +   r   ^^^  ;^^-^^  q(t)f(t,s)dt,  (1.6) 


which  can  be  obtained  by  the  method  of  variation  of  parameters.   With  the 
help  of  these  equations  the  following  lemmas  are  proven: 
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Lenma  1.1  For  each  x  ^  0,   cp(x^s)  is  an  entire  function  of  s  for  vhlch 
the  estimate 

|(p(x,s)l  g  K-^ (1.7) 

1  +  [six 

holds.  Moreover^  (p(x,s)  la  an  even  function  of  s  for  real  s. 

Lemma  1.2  For  each  x  g  0,  f(x;,s)  is  analytic  vlth  respect  to  s  in  the 
upper  halfplane  t  >  0  and  continuoais  dovn  to  the  real  axis.  Moreover,  the 
inequality 

|f(x,s)|  g  Ke"^"",   T  S  0  (1.8) 

holds. 

Lemma  1.5     f(x,s)    satisfies   the   folloving  Inequalities! 

00 

.  -TX  p 

|f(x,s)    -   e^^''!    i  K  ^ /       |q(t)|dt,  t  S   0,  (l.9) 

Is  I         ^ 
'     '  x 

00 

|f(x,s)    -   e^^^l    g  Ke"^^      /     t|q(t)|dt,  t  ^   0,  (l.lO) 

x 

00 

|f(x,s)    -  ise^^^^l    g  Ke""^^      T     |q(t)|dt,  t  S   0.  (l-H) 

X 

The  estimate  (l.9)  is  convenient  for  [s|  -»  oo  and  can  be  applied  when 
X  7^  0.   The  estimates  (l.lO)  and  (l.ll)  are  convenient  for  x  ->  oo.   in  addi- 
tion, from  (l.ll)  one  can  conclude  that 

11m   xf'(x,s)  =  0.  (1.12) 

X  -  0 


All  absolute  constants  depending  only  on  C  (which  can  be  different)  will 
be  denoted  by  K. 
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In  fact,,  as  x  -►  0^ 


"J 

X 


|q(t)|dtg    /   t|q(t)ldt+  7^    /   t|q(t)|dt-0. 


Lemma  l.k     The  function  f(x^s)  for  any  x  is  continuously  dlfferen- 
tiable  vlth  respect  to  s  dovn  to  the  line  t  =  0  with  perhaps  the  exception 
of  the  point  s  =  Q.   The  estimate 

lf(x,s)  -  ixe^^^l  g  -^   e-'^^   T  g  0,  •    (l.lj) 

holds  uniformly  in  x. 

Lemma  1 . 5  For  large  | s ] 

,    ^   sin  sx     /e'  '  \  /-,  -,  i,  \ 

cp(x,s)  =  — -  +  o ,  (l.li)-) 

V  |s|  / 

f(x,s)  =  e^^""  +  o(e"^''),    T  g  0,  (1.15) 

vmiformly  for  all  0  g  x. 

When  s  is  real,  it  is  not  difficult  to  establish  the  relationship  be- 
tween cp(x,s)  and  f(x,s).  Without  further  mention,  ve  will  write  k  instead 
of  s  vhen  s  is  real.   The  solutions  f(x,k)  and  f(x, -k)  =  f(x,k)  for  k  ^  0 
are  linearly  independent  solutions  of  equation  (l.l).   In  fact,  their 
Wronskian 

[f(x,k);f(x,-k)]  =  f'(x,k)f(x,-k)-f(x,k)f'(x,-k)  =  21k     (1.16) 
is  not  zero. 
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In  consequence  of  the  realness  of  cp(x^k) 

(p(x,k)  =  ^  ^f(x,k)M(:¥y  -  f(x,-k)M(k)),  (1.17) 

vhere  M(k)  can  be  determined  with  the  help  of  the  Wronskian  (see  (1.12)): 

M(k)  =  rq3(x,k);  f(x,k)~|  =  11m   ^(x.k);  f(x,k)'l  =  f(0,k)   (l-l8) 

From  this  and  lemma  1.2,  we  conclude  that  M(k)  is  the  limiting  value  of  the 
function  M(s)  =  f(0,s),  analytic  In  the  upper  halfplane  and  such  that 
M(k)  =  M(-k) .  Let  us  introduce  the  notations: 

A(k)  =  lM(k)l,   Ti(k)  =  argM(k),  (1.19) 

A(k)  =  A(-k),   ii(k)  =  -Ti(-k).  (1.20) 

Then  from  lemma  1.^,  we  infer  that  for  large  x 


q)(x,k)  =  ^^   sln(kx  -  Ti(k))  +  o(l),  (l.2l) 


A(k) 
k 

cp'(x,k)  =  A(k)cos(kx  -  ii(k))  +  o(l).  (1.22) 


It  is   therefore  natural  to   call  A(k)   the  limiting  aiiq)lltude  and  Ti(k)    the 
liTirT ting  phase. 

Let  T   >  0.      From  equation   (l.l)    for  f(x, s)    and  the  equation 

-  f"(x,s)    +  q(x)f(x,s)    =   2sf(x,s)   +  s^f(x,s)  (1.25) 

for  f(x,s)  =  df(x,s)/ds,  it  is  not  difficult  to  obtain  the  following  Iden- 
titieai 

00 

f'(0,s)f(0,s)  -  f(0,s)f'(0,s)  =  hla-z      /   |f(t,s)I^dt,      (1.2^) 


16  - 


f'(0,s)f(0,s)  -  f(0,s)f' (0,s)  =  2s   /   f^(t,s)dt.  (1.25) 


From  the  first  one,  we  conclude  that  M(s)  can  vanish  only  for  ct  =  0  or 
T  =  0.   The  second  possihility,  however,  is  excliided  by  the  fact  that  if 
M(k)  =  0  on  the  real  axis_,  (I.I7)  would  imply  that  (p(x,k)  =  0^  and  this 
is  impossible. 

In  the  following,  it  will  be  assumed  that  M(0)  f-   0.   The  vanishing  of 
M(s)  for  s  =  0  is  equivalent  to  the  solution  of  -y"  +  q(x)y  =  0,  y(0)  =  0 
being  bounded  as  x  -►  «>  and  this  happens  only  in  exceptional  situations. 
A  treatment  of  the  case  M(0)  =  0  does  not  present  any  essential  difficulties 
but  only  complicates  the  formulation  and  proof  of  theorems. 

There  still  remains  the  possibility  that  M(s)  =  0  for  a  =  0  and  t  >  0. 
From  the  estimate 

M(s)  =  1  +  0(1)  (1.26) 

for  large  [s],  which  follows  from  formula  (l.l^ti),  we  conclude  that  M(s)  can 
only  have  a  finite  number  of  zeros  of  the  form  s  =  i/c   (n  =  1,  ...,m)  on  the 
imaginary  axis.  When  s  =  s  ,  the  solutions  cp(x, s  )  and  f(x, s  )  satisfy  the 
same  boundary  condition  at  x  =  0  and  therefore  they  are  proportional 

f(x,sj  =  f'(0,s^)cp(x,s^).  (1.27) 


From  this  and  (I.25),  it  follows  that 


/    Ux,s^)J"dx=-^,.y,3  )  ,  (1.28) 

a    I-      -I         n      n 
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and  this  implies^,  in  particular^  that  all  the  zeros  of  M(s)  are  simple. 

The  above  results  can  be  sunMarized  in  the  following  form. 

Lemma  1.6.   The  function  M(s)  Is  analytic  in  the  upper  halfplane  and 

has  there  a  finite  number  of  simple  zeros  s  =  ±k   ,    K     >  0,  (n  =  1, ...,m). 
^ n     n   n  7    ^  / 

For  large  [s|,  the  estimate  (1.26)  holds.  The  function  M(s)  is  continuous 
down  to  the  real  axis  vlth  the  possible  exception  of  the  point  s  =  0.  Fur- 
thermore^ sM(s)  is  continuous  everywhere.   This  last  assertion  follows  from 
lemma  1.^. 

2.   Expansion  Theorem 

The  differential  equation  (l.l)  together  with  the  boundary  condition 
defines  a  self-adjoint  operator  In  o^p(0^«>).   This  operator  can  be  gotten 
by  the  extension  of  the  symmetric  operator^  defined  by  (l.l)^  acting  on  the 
twice- continuously  differentiable  functions  satisfying  the  boijndary  condi- 
tion and  vanishing  identically  outside  some  finite  interval.  We  will  denote 
this  operator  by  L. 

Consider  the  kernel 

^  f(y,  A) 
^  M(  A) 


Rj^(x,y)  =  R^(y^x), 


0  ^  arg  /X  ^  It, 


(2.1) 


which  is  defined  for  all  complex  X  with  the  exception  of  a  finite  number 
of  points  on  the  negative  real  axis  corresponding  to  the  zeros  of  M(  ^Jx) . 
By  virtue  of  (I.18);,  it  is  not  difficult  to  verify  that  the  kernel  R  (x^y) 

A 


18 


is  a  solution  of  the  equation 


f-  ^  +  q(x)j  R^(x,y)  -  XR^(x,y)  =  5(x  -  y) ,  (2.2) 

and  also  satisfies  the  boundary  conditions: 

\(0,.y)  =  R^(x,0)  =  0.  (2.3) 

In  consequence  of  (l.?)  and  (1.8)^  we  have 


K  (x,y)  I  ^  K  —^ e"^l''"^l,T  =  Im  yi  >  0  (2.1+) 

1  +  I  /x\x 

for  complex  X  and,  hence,  the  kernel  IL  (x,y)  determines  In  a(f_(0,i»)  a  bounded 
operator,  namely,  the  resolvent  operator 

Rj^  =  (L  -  XI)"^.  (2.5) 

The  singularities  of  the  resolvent  R  In  the  complex  X-plane  consist  of 

a  cut  along  the  positive  real  axis  and  a  finite  number  of  simple  poles 

2 

X  =  -K      (n  =  l,...,m)  on  the  negative  real  axis.   The  continuous  and  dis- 
n     n      f        }    /  o 

Crete  portions  of  the  spectrum  correspond  to  the  cut  and  poles,  respectively. 
The  jump  in  the  resolvent  across  the  cut  and  the  residues  at  the  poles  deter- 
mine the  spectral  function  of  the  operator  L.  We  come  now  to  the  following 
completeness  theorem  for  the  elgenfunctlons  of  the  operator  L. 

Theorem  2.1.   The  functions  cp(s,k)  (k  &  O)  and  cp  (x)  =  qp(x,iK  )  form 

j^  n 

a  complete  orthogonal  set  of  functions.   The  completeness  relationship  is 
given  by 
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Z  ^n%^-K^^)   ^   f  f    ^(^^^)  M(k)M(-k)  <P(y'^)>^^^  =  ^(^-y)        (2.6) 
n=l  "i 


in  which  C  =  2±K  f'{0,lK   )/M(i/c  )   (see  (1.28)). 

Formula  (2.6)  can  be  deduced  without  any  reference  to  the  general  aspects 
of  operator  theory.   Let  f(x)  be  a  twice  continuously  differentiable  function 
vanishing  for  large  x  and  in  the  neighborhood  of  x  =  0.  Then 

g(x)  =  -f"(x)  +  q(x)f(x)  (2.7) 

is  contlnuoiis  and  vanishes  identically  outside  some  finite  Interval  not  con- 
taining the  origin.  From  (2.2)  and  (2.7),  it  follows  that 

J    R^(x,y)f(y)dy  =  -  i  f(x)  +^   J    R^(x,y)g(y)d7.  (2.8) 

o  o 

If  we  take  the  integral  of  both  sides  of  (2.8)  along  a  large  circle  [x]  =  N, 
the  contribution  from  the  second  term  on  the  right  side  of  (2.8)  as  N  "^  <» 
will  be  zero  by  virtue  of  (l.l^),  (1.15)  and  (1.26).   Th\is  we  get 

00 

lim   ^    (S         r  r  R^(x,y)f(y)dy  dX=-f(x).         (2.9) 
^^~        xilNl    o 

On  the  other  hand,  if  we  take  the  integral  of  the  left  side  of  (2.8)  along 

a  path  7  consisting  of  a  curve  encompassing  the  cut  on  the  i^al  axis  and  a  large 

circle  Ixl  =  N,  we  obtain 
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il  f  if  \(x.y)^(y)^y  Ux  ^  ^    f     If  R^(x,y)f(y)dy  I  ^  + 


+  2^    /    ^'{    f     h^lO^^'^^    -   ViO^^^y)^^^)^^ 


o 


m 


=    Zl    ^^^    {   f   R.(x.y)f(y)dy 


1  u  ^ 

n=l  I  ^^ 


n 


By  taking  into  account   (2.1)   and  (I.I7)    and  passing  to  the  limit   as  N -♦  <»   _, 
ve  find  on  the  basis   of  (2.9) 


00  {     °° 

o  1^0 


j-^  q)(y.k)f.(y)dy   >    + 


m  00 


+    X!      r    C  cp^(x)cpjy)f(y)<iy 
n=l     '^ 


Finally,  by  virtue  of  the  fact  that  the  functions  f(x)  are  dense  inXp(0,oo), 
formula  (2.6)  is  obtained. 

The  functions  ilr^'*"'' (x^k)  =  cp(x,k)/M(k)  and  i|f  (x)  =  /c~  cp(x)  form  an 
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orthonormal  system.  Hovever,  the  functions  ^^  (x^k)  are  not  square-inte- 
grable  and,  hence,  are  not  elements  of  Hilbert  space  and  eigenfunctions  in 
the  usual  sense.  To  attach  meaning  to  them  vhlle  still  remaining  in  the 
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framework  of  Hllbert  space,  one  may  treat  them  as  kernels  of  transforma- 
tions, vhich  diagonallze  the  operator  L.   Thus,  the  transformation 

00 

T'"^^g  =  G:   G(k)  =   /  g(x)t'^'^^(x,k)dx  (2.10) 

o 

carries  any  function  g(x)  in  «!f  p(0,oo)  into  a  fimction  G(k)  for  which 

00 

|G(k)|  k;  dk  <  ».  Moreover,  if  Lg(x)  belongs  to  ^p(0,<»),  then  Lg(x)  goes 

00 

into  k  G(k)  and  the  integral   /   |G(k)k  |  K  dk  converges.   In  the  following^ 

o 
the  space  of  functions  G(k)  with  the  scalar  product 

00 

(G,G^)  =  I  /  G(k)G^(k)k^dk  (2.11) 

o 

will  be  denoted  by    iC.}    and  the  space  of  square-integrable  f\mctions  g(x), 

which  was  previously  called    jig    (O.oo),   will  be  denoted  by  ^    .      The  trans- 

formation  T         acts  from   ^     into    X^.      The   adjoint  transformation  T^ 

acts   from   iZ,    into   aC      according  to  the  formula 
k  X 


00 


g  =  T^'*'^*G:      g(x)   =  -     /'     G(k)ilf^'^nx,k)k^dk.  (2.12) 


The  orthogonality  of  the  eigenfunctions  t^^      (x,k)  can  thus  be  expressed  by 
\ising  the  new  terminology  in  the  following  way: 

T^-^^T^-^)"  =  I^  (2.13) 
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Here,  I  denotes  the  identity  operator  in  ;^  .   Relation  (2.6)  can  nov  be 
written  down  in  the  following  form: 

Here,  I  is  the  projection  operator  onto  the  proper  subspace  of  the  opera- 
tor L  corresponding  to  the  continuous  spectrum.   The  superscript  x  indicates 
that  this  operator  acts  on  Jf  .  A  transformation  of  the  type  T    can  be 
associated  with  aay  solution  of  equation  (l.l)  satisfying  a  zero  boundary 
condition.   Let  X(x,k)  be  solution  of  (l.l)  such  that  X(0,k)  =  0  and 
X'(0,k)  7^  0  for  all  k.   Consider  the  transformation 


Tg 


=  G!   G(k)  =   /  g(x)x(x,k)dx.  (2.15) 


The  transformations   T  and  T         are  related  by  the   foiroula: 

T  =  N^T^"^)  (a.  16) 

where  IP  is  a  'normalizing  factor';  N  is  an  operator  in  ^  ,  which  multi- 
plies  by  the  function  N(k)  =  X'(0,k)/i|f    (0,k).   The  completeness  relation 
and  orthogonality  condition  can  be  written  down  in  terms  of  T  in  the  fol- 
lowing fashion: 

T^W^  =  I^  ,   TT*W^  =  l\  (2.17) 

where 

W^  =  (/)"^(N^*)"^.  (2.18) 

The  somewhat  formal  arguments  at  the  end  of  this  section  will  turn  out  to 
be  useful  in  Sec.  5« 


-  25  - 


5.   Aaymptotic  behavior  of  the  solution  of  the  Shroedlnger  equation  for 
large  time. 

The  expansion  theorem  for  the  operator  L  deduced  In  the  preceding  sec- 
tion enables  one  to  solve  by  Fourier  methods  time-dependent  equations  In- 
volving the  operator  L.  We  will,  in  particular,  be  interested  in  the  be- 
havior of  solutions  of  the  Schroedlnger  equation 

i  ^4f^  =  Lf(x.t),  f(x,t)  1^^^  =  f Jx)  (3.1) 

for  large  | t | . 

By  expanding  f  (x)  vlth  respect  to  the  eigenfunctions  of  L,    one  can 
represent  the  solution  of  (5-1)  in  the  following  way: 

f(x,t)  =  -   r  F(k)  tC^^   e"^^  Vdk  +     ^     F     TUJ   e  "^  . 
-  o 

(5v2) 

Here  i|f(x,k),  -((r  (x)  form  an  arbitrary  set  of  orthonormal  eigenfunctions  of 
L     i|f(x,k)  and  ijr  (x)  may  differ  from  the  functions  f^      (x,k)  and  f   (x) 
dlsciissed  in  Sec.  2  by  a  factor  of  modulus  one   and 

00  oe 

F(k)  =  J   f^(x)t(x,k)dx,  F^  =  jT  f^(x)tjx)dx.  (5.5) 

o  o 

In  particular,  one  may  take  as  the  set  \)f(x,k)  the  functions  \)/^   (x,k)  or 


t^'^(x,k)  =  i|r^''"\x,k)M(k)/M(-k)  =  t   (x,k) .   The  behavior  of  f(x,t)  for 
large  |tl  can  be  analyzed  on  the  basis  of  the  following  lemma: 


-   2k 


Lemma  ^-l*      Let  F(k)  be  an  arljitrary  function  la  ^    ,    I.e. 

00 


o 


and  set 


f^*^(x,t)   =      r    F(k)   t^*^(x,k)   e"^^  ^  k^dk,  (3.5) 


2 


Then 


g(x.t)  =    J    F(k)  ^ii^  e-^^  ^  k^dk.  (5.6) 

o 

00 

lim         f    |f^*^(x,t)    -  g(x,t)|^dx  =  0.  (3.7) 

->    +  00  '-' 


t-±"    ^ 


It   is  sufficient  to  prove  this  assertion  for  a  set  of  functions  F(k) 
vhich  is   dense  in  )C ,)    for  an  arbitrary  function,    the  theorem  then  follovs 
by  completion.      We  will  assume  that  F(k)    is  differentlable   and  does  not 
vanish  in  the  interval  0<a^kg    p<oo.      The  inteirvals  of  integra- 
tion in  (3-5)   an<3.   (3-6)   are  then  finite   and  do  not   contain  the  point 
k  =  0.      For  definiteness,   we  will  assxjme  that  t  ^  »    .      By  virtue  of  the  fact 
that  i|;^      (x,k)   and  (sin  kx)yk  are  iiniformly  bounded  with  respect  to  all  x 
in  0  <  X  <  <»  and  all  k  in    [a,  p] ,    the   functions   f^      (x^t)    and  g(x,t)   tend  to 
zero  as  t  -»•  00     uniformly  in  x.      Thus  the  integral  in   (3-7)   vanishes  as 
t  -»•  00     for  any  finite  Interval  of  integration.     We  still  have  to  show  that 

00 

/    |f^   -  g|  dx  converges  for  arbitrary  A. 
A 

From  (1.17)  and  the  definition  of  t   (x,k)^ 
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T^'W)  -^%^ 


1   -ikx 
e 


M(k 


2ik      [M(-k 


-1 


+  R(x,k),       (3.8) 


whei^  for  x  >  0 


|R(x,k)|  S  K  4   /   |q(t)|dt 
k  J 

X 


(5.9) 


by  virtue  of  (1.9)  sud   (1.8).   In  consequence  of  this  estimate, 


00      oo 


2 

F(k)R(x,k)e"^^  Vdk 


dx  g 


A    a 


SK   /   d^  (  /   k(t)|dt  J    /   |F(k)|\' 
A     ^  X         ^   a 


^dk  /   ^. 

^   k 
a 


g  K'  r  iq(t)idt  r  t|(i(t)idt 


(5.10) 


and  for  sui'ficiently  large  A  the  integral  containing  R(x,k)  can  "be  made  as 
small  as  desired  uniformly  in  t.  We  still  mujst  look  at  the  behavior  of  the 
integral  of  the  basic  term  of  (5.8) 


^  /a.^     r  ^     r    nr^  \    -ikx  -ik  t^, 
Q  (t)  =   /   dx    /   G(k)e    e     dk 


(5.11) 


a 


as  t  ->■  CO 


Ffk)  Mfk 
The  function  G(k)  =  -gr-^  rTr- 


^ 


-  1 


k  by  virtue  of  lemma  1.6  is 


.B, 


finite  and  continuously  differentiable.  Now  Q.(t)  =  lim  Q.(t)  where 

B   -*-   00 
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<(t) 


dx 


'G(k)e-i^e-^^  *dk 


a 


P    P 


^-I(k2_e2)t  ,-l(k-l)A_^-l(k-l)BJ  ^  ^     _ 


dk  /  dt   ^G(k)GriTe  ^^-  -^  ^^  ^ ^^.^^^^^ ^  =  J^(t)-J^(t), 


a    a 


Here 


P    P 


-i(k^-l^)t  -l(k-t)B 
J^(t)  -J^J   d^G(k)G(Ty^ ^^-^^ , 

a    a 

and  because  of  the  singularity  in  the  denominator,,  the  inner  integral  is 
taken  els  a  principal  value,  well-defined  in  consequence  of  the  differentia- 
bility of  G(k). 

We  transform  J^{'^)   to  the  following  form.: 


Jg(t)  =  /  dkG(k) 


I-  B  2  2 

r        GnTe-^^^-^^^-GTky  ^-i(k-e)B 


a 


i_a 


P  ^-i(k-^)B 
+  ^^j   i(k-i)  ^' 


(5.12) 


The  first  term  in  the  integrand  is  continuous  for  k  =  I  and  vanishes  as 

B  -*  00  by  the  Riemann-Lebesgue  theorem.   The  integral  of  the  second  term 

tends  to  jt  as  B  -►  »  and  therefore,  we  find  that   lim   J.o(t)  does  not 

3  B-.C.   B 

depend  on  t  and  is  equal  to  «   /   |G(k)|  dk. 

a 
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J  (t)  can  also  be  represented  as  the  sum  of  two  terms  one  of  which. 

vanishes  as  t  -►  »  uniformly  In  A.   The  other  term  has  a  finite  limit 

P 
Independent  of  A  as  t  -►  «  equal  to  rt   /   lG(k)|  dk..   From  this  it  follows 

a 
that  Q. (t)  -*  0  and  the  lemma  is  therefore  proven. 

Let  us  now  return  to  the  investigation  of  the  behavior  of  the  solu- 
tion f(x^t)  of  the  Schroedlnger  equation.   If  f  (x)  is  orthogonal  to  the 
eigenfunctions  of  the  discrete  spectrum  of  L,    the  sum  in  the  second  term  of 
(5-2)  does  not  appear.   Let 

00 

F^*^(k)  =  r    f^(x))|/^*\x,k)dx.  (5.15) 

o 

These  are  the  fvinctions  which  occur  in  formula  (5-2)  for  f(x,t),  the  \(f^~'^(xjk) 
having  been  selected  for  the  functions  ilf(x^k) .   It  is  evi  .ent  that 

00 

(±)/   ^^    2   r  ^(±)/,  N  sin  kx  -ik^t  .  2^,  /'^^k^ 

g^   (x,t)  =  -  /   F^  '(k) —  e      k  dk  (5.1^j 


is  a  solution  of  the  Schroedlnger  equation  with  an  operator  L  ,    associated 

o 
with  the  equation  L  y  =  -y"  =  k  y  and  the  boundary  condition  y(0)  =  0^  i.e. 

lMi!^.L„g(*)(x,t),  (5.15) 

where  „ 

g(*'(x,o)  .  g(*'(x)  .  f  /  /*-h^}  iii^  k^dJ..      (5.16) 

o 

If,  in  analogy  with  Sec.  2,   we  Introduce  a  unitary  transformation  T  of  ;C 
Into  X]5.  given  by 
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T^g  =  G:   G(k)  =  J    g(x)  ^i£^  dx,  (5.IT) 

o 

then  g^~  (x)  -will  be  expressed  in  terms  of  f  (x)  by  the  formula: 

g^-^(x)  =  W-)f  fx).  (3.18) 


o  '      o     o' 


From  lemma  5-1  follovs 

Theorem  '^.1.      If  f  (x)  Is  orthogonal  to  the  elgenfunetlons  of  the  dis- 
crete spectrum  of  L^  then  the  solution  of  the  Schroedlnger  equation  (3.I) 
as  t  -*  ±  00  behaves  like  the  solution  of  the  Schroedlnger  equation  (3«1^)  vlth 
initial  data  g   (x)  given  by  (3'l8),  in  the  sense  that 

/   |f(x,t)  -  g^-^(x,t) l^dx  ^  0  as  t  -  ±  "  .  (3.19) 


k-.        Transformation  operators. 

In  the  following  ve  will  need  the  representation  of  a  solution  of  (l.l) 
with  a  potential  q(x)  in  terms  of  solutions  of  the  equation  with  other  poten- 
tials and,  in  particular,  with  q(x)  =0,  i.e.  in  terms  of  trigonometric  func- 
tions. 

The  silirplest  expression  of  this  kind  can  be  deduced  in  the  following 

manner.   On  the  basis  of  (1.9)  for  x  ^  0,  the  function  h(x,s)  =  f(x,s)-e 

is  square-integrable  in  s  along  any  line  parallel  to  the  real  axis,  and  in 

00 

|h(x,  a   +  It) I  da  =  0(e    ).   By  a  theorem  of 
-  00 

Tltchmarsh, 
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A(x,y)   =^   f      (f(x,k)-e^^)  e'^^dk  =  0,    x  >  y.  (4.1) 

-    00 

Inverting  this  Fourier  transform^    one  thus  obtains   for  f(x, s) 

00 

f(x,s)   =   e^^^  +    J    A(x,y)e^^^dy,    x  S  0,  (4.2) 

X 

vhere  A(x,y)  Is  square -Integrable  in  y  for  x  ^  0.     With  certain  modifications^ 
this  procedure  yields  an  expression  for  cp(x, s)  valid  for  all  s 

X 

cp(x,s)  =  ^^  +  J    K(x,y)  ^iS_M  ay.  (4.5) 


o 


However,  the  derivation  gives  alcaost  no  information  concerning  the  kernels 
K(x_,y)  and  A(x,y).   Equivalent  equations  for  these  kernels  can  be  deduced  by 
substituting  expressions  (4.2)  and  (4.5)  for  f(x^s)  and  9(x,s)  into  (1.6)  and 
(1.5),  respectively,  and  by  eliminating  the  trigonometric  functions.   Doing 
this,  one  obtains  the  equations 

x+y  x+y      x-y 

~2~  2        2 

K(x,y)  =  I   r  q(t)dt  +   r   dt    r  dz  q(t+z)K(t+z,t-z),   x  g  y,    (4.4) 

x-y  x-y      o 

2  2 


y-x 
00         00     2 

^(x^y)  =  i   r  q(t)dt  -   I        dt       I     dz  q(t-z)A(t-z,t+z),  y  ^  x.    (4.5) 

x+y        x+Z     ° 

2  2 


By  then  solving  these  equations  by  the  method  of  successive  approximations. 
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ve   get  the   inequalities 


x+y 
2 


|K(x,y)|  §  I      J     |(i(t)ldtexp        Jt\ci(t)\dt         ,  (4.6) 


x-y  o 

2 


A(x,y)|gi        r     |q(t)|dtexp        rtlq(t)|dt  {k.l) 


2 

x+y 
2 


for  K(x,y)  and  A(x,y).   From  these  integral -equations^  it  also  follows  that 
K(x,y)  and  k{x,j)    are  differentiable  and  estimates  for  their  derivatives  can 
be  determined.   For  example^ 

00  00 

|^A(x,y)  +  Jq(^)|  ^K  J     |q(t)|dt  J     k(t)ldt  .      (1^.8) 

X  x+y 

2 

A  similar  inequality  holds  for  dA(x^y)/Sy. 

By  virtue  of  (4.6)  and  {^.l),    the  theory  of  Volterra-type  integral 
equations  may  be  applied  to  integral  equations  having  K(xjy)  and  A(x^y)  as 
kernels.   Thus^  if  g(x)  belongs  to  a  given  class  of  functions  (which  we  do 
not  specify),  then  the  equations 

X 

Ugf(x)  =  f(x)  +  J    K(x,y)f(y)dy  =  g(x),  (4.9) 

o 

00 

Vgf(x)  =  f(x)  +  J    A(x,y)f(y)dy  =  g(x)  (4.10) 


51 


have  solutions  vhich  can  be  represented  in  the  form: 


f(x)  =  g(x)  +  J   K(x,y)g(y)dyH  V^^ix) ,  {^.H) 

o 

00 

f(x)  =  g(x)  +  J    A(x,y)g(y)dy=  Vg^g(x)  .  (^^.12) 


In  connection  with  this^  the  kernels  K(x,y)  and  A(x,y)  will  have  estimates 

similar  to  (^.6)  and  (k.j).     A  more  precise  definition  of  the  operators  U 

and  V  as  operators  in  Hilbert  space  will  be  given  in  the  subsequent  sections. 
B 

Relationships  analagous  to  (^.5)  exist  between  any  two  solutions  of  (l.l) 
with  different  potentials.  Let  cp^(x,k)  and  cpp(x,k)  be  solutions  of  (l.l)  with 
potentials  q, (x)  and  q  (x),  respectively:  two  relations  of  type  (^.3)  can  then 
be  derived: 

/  ,  V   tt(1)  sin  kx       c  ,  \   tt(2)  sin  kx  /,  ,^\ 

(p^(x,k)  =  U^   ^—  .   cpgCxA)  =  U^   ^—  .  (^.15) 

Inverting  the  second  of  these  and  substituting  the  result  in  the  first,  we  get 
OT,    expli  citly 

X 

(p^(x,k)  =  q>2(x.k)  +  /  K(x,y)cp2(y/k)dy  .  (^-15) 


It  is  also  possible  to  derive  an  inequality  for  K(2,y)  similar  to  (k.6). 
Finally,  one  may  also  relate  two  different  solutions  f  (x,k)  and  f  (x,k)  but 
such  an  expression  will  not  be  reqxiired  in  the  following. 
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5-   General  theoiy  of  transformation  operators. 

We  will  consider  the  operator  L  Introduced  in  Sec.  2  as  one  of  the  func- 
tional representations  of  an  abstract  operator^  vhlch  will  also  be  denoted 
by  L.   The  operator  L  of  Sec.  2  will  now  be  denoted  by  L  and  the  space  pf   In 
which  it  operates  will  be  called  the  coordinate  representation  or  x-represen- 
tatlon.  Another  representation  to  be  considered  is  the  so-called  momentum  or 

k-reprasentation.   This  space  is  of  type  X^y   and  it  is  defined  by  the  condl- 

2 
tion  that  in  it  L  is  an  operator  which  multiplies  an  element  by  k  : 

L^F(k)  =  k^(k).  (5.1) 

Both  representations  are  related  to  one  another  by  the  unitary  transformation 
T  introduced  In  Sec.  5-  In  other  words^  to  each  element  f(x)  e  X,      there 
corresponds  an  element  F(k)  belonging  to  JC.  ! 

00 

F(k)  =  T^f(x)  =  ^    f(x)  ^^£^  dx  (5.2) 

o 

V  k      X 

and  each  operator  A  in  ;i^  is  converted  into  an  operator  A  =  T  A  T*.  For 

example^  in  the  k-representation^  L  is  given  by 


L^(k)  =  k^(k)  +  I  r  v(k,^)F(^)^^d^,  (5.5) 


vhere 


/,     .\  r    sin  kx     /    %    sin  tx    ,  /c   i,  ^ 

V(k,t)  =     /     — Y~  'i(^)  — T~  ^'  ^^'    > 


o 


The  passage  from  the  momentum  back  to  the  coordinate  representation  is  effected 
by  using  the  transformation  T*^  so^  for  example,  the  operator  vhlch  multiplies 


33 


by  a  decreasing  function  n(k)  goes  into  an  Integral  operator  in  the  x-repre- 
sentatlon  vlth  a  kernel 

00 

„/    \   2   r  sin  kx  ^,.  ^    sin  ky  ,  2,,  fc:   ^\ 

n(x,y)  =  -  /  i^  fXk)  j~  k  dk.  (5.5) 

o 

In  Sec.  2,  the  elgenfunction  expansion  theorem  for  L  was  written  down  in  terms 
of  the  operator  T  and  involved  various  spaces.   This  is  inconvenient.   In  place 
of  T,  we  will  utilize  operators  acting  in  one  and  the  same  space  and  which 
therefore  can  be  given  by  one  of  their  representations.   These  operators 
are  defined  in  the  coordinate  and  momentum  representations  in  the  following 
way: 

IT"  =  T  T*  in  the  k-representation 

,   °  I  (5-6) 

U^  =  T^T   in  the  x-representation  J 

Such  operators  will  be  called  transformation  operators.   In  all  subsecjuent 
discussions,  they  replace  completely  the  eigenfunctions  of  the  continuous 
spectrum.   Instead  of  a  differential  equation  for  the  eigenfunctions  X(x,k), 
one  has  for  the  operator  U,  determined  by  them,  the  equation 

LU  =  UL^.  (5.7) 

The  completeness  condition  and  orthogonality  of  the  eigenfunctions  In  terms 
of  U  are  expressed  by 

UWU*  =  I  ,       U*UW  =  I  (5.8) 

(see  (2.17)).   The  operator  W  which  appears  in  these  formulae  is  an  operator 
which  multiplies  by  the  function  W(k)  in  the  momentum  representation.   The 
convenience  of  (5.7)  snd   (5.8)  is  in  not  having  to  write  indices  indicating 
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the  spaces  In  which  the  operators  act  -  all  operators  are  in  one  and  the 
same  representation.   In  the  many  problems  related  to  the  operator  L^  the 
use  of  different  transformation  operators  is  advantageous.   Thus^  the  opera- 
tors U^"  obtained  from  the  transformations  t    by  means  of  (5-6)  turn  out 
to  be  useful  in  studying  the  asymptotic  behavior  of  the  solution  of  the 
Schroedinger  equation 

i^=Lz(t),   z(t)  ,_  =  z^  (5.9) 


"3t" 


t=o    o 


for  large  |t|.   The  solution  can  be  written  in  the  following  form: 

z(t)  =  e-^^^z^,  (5.10) 

and  theorem  3.1  can  now  be  reformulated  as 

Theorem  5.1.   If  z  is  orthogonal  to  the  eigenfunctions  of  the  point 
spectrum  of  L,  then  for  the  solution  of  the  Schroedinger  eguatlon  the  limiting 

conditions 

IL  t 
11m   e  °  z(t)  =  z^  (5.11) 

t  -»  ±00  ~ 

hold  where 

z  =  u'^''"^*z  ,    z  =  u''"''*z  .  (5.12) 

If  we  denote  the  'normalizing  factor'  2?elatlng  U^'*''^  and  U^""^  by  Si 

u(-)  =  u^+)s,  (5.15) 

then  z  and  z  in  formula  (5.12)  are  related  to  one  another  by: 

z^  =  Sz_.  .  (5.1^) 

It  is  not  difficult  to  show  from  the  definition  of  u'''*"'^  and  (5-15)  that  in 

the  momentum  representation  this  is  an  operator  which  multiplies  by  the  function 
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S(.)  =  ^  .  (5.15) 

S(k)  obviously  has  an  absolute  value  equal  to  one.   From  theorem  '^.1,    it 

follows  that  there  exists  a  solution  of  the  Schroedinger  equation  which  is 

-iL  t 
asymptotic  to  the  vector  z  (t)  =e     z  ast'-^-oo^z  being  an  arbitrary 

element.   In  addition,,  this  solution  for  t  =  0  will  be  orthogonal  to  the 

eigenfunctions  of  the  point  spectrum  of  L  and  consequently  will  behave  like 

-iL  t 
e     z  ast-*-oo^z  being  constructed  from  (5«1^)  using  the  operator  S. 

Thus^  the  general  formulation  of  the  scattering  problem  described  in  the 
introduction  turns  out  to  be  valid  for  L,    and  S  plays  the  role  of  the  scat- 
tering operator.  We  have  also  shown  that  S  is  unitaiy  and  commutes  with 
the  energy  operator. 

The  fundamental  aim  of  the  survey  is  to  establish  the  relationship 
between  L  and  S.   The  operator  U^  =  I  +  K,  defined  in  the  coordinate  repre- 
sentation  by  (^.9)^  will  play  an  important  role  in  this.   By  the  definition 
of  this  section^  this  operator  is  a  transformation  operator.   In  fact^  if  we 
assume  the  transformation  operator  to  be  an  integral  operator^  then  its  ker- 
nel in  the  x-representation  is  obtained  by  expanding  the  corresponding  solu- 
tion X(x,k)  with  respect  to  .(sin  kx)/k: 

to 

U(x.y)  =  f  f   W^  ^^  k^dk  = 
o 

00 

=  f  f    (X(x.k)  -  ^^)  ^^  k^dic  +  6(x-y).      (5.16) 
o 

In  exactly  this  way,  the  kernel  of  the  operator  U  -  I  with  a  corresponding 
solution  cp(x,k)  was  deduced  in  Sec.  k-.      One  can  give  an  exact  meaning  to  this 
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crude  argument.  In  consequence  of  the  fact  that  the  function  M(k)  is  the 

normalizing  factor  for  the  solution  cp(x,k),  with  the  help  of  which  U  is 

n 

constructed,  the  operator  W  appearing  in  the  formula  for  U  of  type  (5.8) 

B 

is  given  hy  the  function 


^^^)  =  M(k)M(-k) 


(5.17) 


A  characteristic  feature  of  U^  is  that  it  is  a  Volterra  operator  in 
the  coordinate  representation.   This  property  is  closely  related  to  the 
fact  that  the  potential  q(x)  is  diagonal  in  this  representation  (i.e.  it  is 
a  multiplication  operator).   In  fact,  the  triangularity  of  the  kernel  K(x,y) 
is  a  consequence  of  the  fact  that  cp(x, s)  is  an  entire  function  of  s  which  is 
due  to  the  diagonal  property  of  the  potential.   Conversely,  let  there  exist 
for  some  operator  L  =  L  +  V  a  transformation  operator  U  of  Volterra-type 
in  the  x-representation. 


Ug(x,y)  =  &(x-y)  +  Ti(x-y)K(x,y) 


(5.18) 


where 


Tl(t)  = 


1  for  t  >  0 


0  for  t  <  0. 


(5.19) 


From  (5.7),  we  find 


V(I  +  K)  =  KL  -  L  K. 
^     '00 


(5.20) 


The  operator  B  =  KL  -  L  K  has  the  kernel 


B(x,y)  =  26(x-y)  —  K(x,x)  +  Ti(x-y) 


dx 


^^(x,y)  _  a^(x,y) 


5x         Sy 
Multiplying  (5. 20)  by  (l  +  K)~  ,  we  deduce  by  virtue  of  (5.21)  that 


(5.21) 


V(x,y)  =  6(x-y)2  -~~  K(x,x)  +  Ti(x-y)C(x,y) 


dx 


(5.22) 
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and  since  V  has  to  be  self-adjoint^  that 


C(x,y)  =  0,   V(x,y)  =  6(x-y)2  —  K(x,x), 


(5.23) 


i.e.^  L  is  defined  by  an  equation  of  type  (l.l).  A  precise  statement  of  the 

above  conclusions  is  contained  in  the  following  theorem: 

Theorem  5.2.  Let  L  be  an  operator  defined  in  the  x-representation  by 

the  differential  expression  L  y  =  -y"  and  the  condition  y(0)  =  0.   If  there 

exists  a  self-adjoint  transformation  L  =  L  +  V  and  Yolterra  transformation 

U  =  I  +  K,  whose  kernel  is  differentlable  and  if  (5.7)  Is  satisfied,  then 
B 

In  the  x-representatlon_,  V  is  an  operator  which  multiplies  by  the  function 
q(x)  =  2^K(x,x).  (5.2^) 

If  the  kernel  of  the  transformation  is  not  differentlable,  then  it  is 
still  possible  to  define  L  by  the  differential  expression  Ly  =  -y"  +  q(x)yj 
in  this  case,  q(x)  is  a  generalized  function  such  as  the  derivative  of  an 
ordinary  function,  and  in  particular,  may  have  6-function  singularities. 

From  theorem  5.2,  ve  have  simultaneously  shown  that  the  kernel  K(x,y) 
satisfies  the  equation 


q(x)K(x,y)  =  -    p 


5^(x,y)  _  S^(x,y) 


a/ 


X  >  y. 


(5.25) 


In  fact,  by  substituting  (5.23)  into  the  left  side  aad  (5'2l)  into  the  right 
side  of  (5.20),  we  find 

q(x)&(x-y)  +  Ti(x-y)q(x)K(x,y)  =  q(x)6(x-y) 


+  r\{x-j) 


a^(x,y)  _  a^(x,y) 


Bx^ 


V 


(5.26) 
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and  from  this  follovs  equation  (5-25) .   Of  course^  in  the  general  case^  when 
K(x^y)  does  not  have  derivatives^  it  is  a  generalized  solution  of  this  equa- 
tion.  On  the  basis  of  this  equation  for  K(x^y)  and  the  conditions  (5.2^) 
and  K(x,0)  =  0  (this  follows  directly  from  the  definition  of  K(x,y)),  it  is 
not  difficult  to  derive  (Jj-.^),  already  used  earlier. 

6.   The  functions  W(k)  and  S(k),  their  properties  and  relationship. 

In  this  section^  we  will  study  the  functions  W(k)  and  S(k)  which  were 
introduced  in  the  preceding  section  by  the  formulas: 

^^^)  -M(k)i(-k)  (-<^<°°)^  (6.1) 

S(k)=^     (-~<k<»),  (6.2) 

and  we  will  establish  the  relationship  between  them.   In  addition  to  the  pro- 
perties of  M(k)  assembled  in  lemma  1.6^  we  still  require  one  more  property: 
Lemma  6.1.   M(k)  has  the  representation 


M(k)  =  1  +   /  F(t)e^^*dt,  (6.5) 


where  |F(t) |  is  an  Integrable  function. 

The  proof  follows  immediately  from  the  definition  of  M(k)  and  the  expres- 
sion (^.2).   These  imply 

00 

M(k)  =  f(0,k)  =  1  +  rA(0,t)e^^*dt,  (6.4) 

o 

but  lA(0^t)|  is  integrable  because  of  the  inequality  (4,7). 
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From  (6.1)  and  (6.2)  and  by  using  Wiener's  theorem  on  the  Fourier  trans- 
form of  absolutely  Integrahle  functions^  one  can  deduce  the  following  lemmas: 
Lemma  6.2.  W(k)  possesses  the  following  properties: 

1)  W(k.)  is  a  positive  even  function 

W(k)  >  0,   W(-k)  =  W(k)    (-oo<k<=o);  (6.5) 

2)  W(k)    -   1  has   an  absolutely  Integrable  Fourier  transform 

00  00 

W(k)  =1+1     H(t)e^^^dt  =1+2     /     H(t)co3  kt  dt.  (6.6) 

-00  o 

Lemma  6.5  The  function  S(k)  possesses  the  following  properties: 

-1 

1)  |s(k)|  =  S(0)  =  S(<»)  =  1,  S(-k)  =  S(k)  =  Slkj    ;  (6.7) 

2)  S(k)    -   1  has   an  absolutely  integrable   Fourier  transform; 

00 
S(k)   =  1  +      r    F(t)e"^^^   dt;  (6.8) 


5)      arg  S(k) 


-00 

00 


=   -4-irtm,        m  s  0,  (6.9) 


where  m  is  the  number  of  discrete  eigenvalues. 

The  last  property  is  obtained  on  the  basis  of  Cauchy's  theorem  for  the 
number  of  zeros  of  an  analytic  function  by  taking  into  account  that 
arg  S(k)  =  -2  arg  M(k) .   The  properties  of  the  functions  W(k)  and  S(k)  enu- 
merated in  these  lemmas  are  characteristic.  For  if  they  are  fulfilled,  there 
exists  a  \mlque  function  M(s)  which  is  analytic  and  bounded  in  the  upper  half- 
plane_,  and  behaves  asymptotically  like  M(s)  =  1  +  o(l)  for  large  [s|.  Fur- 
thermore, M(s)  has  a  finite  number  of  simple  zeros  s  =  Ik  (n  =  1,  ...,m), 
and  on  the  real  axis,  (6.1)  and  (6.2)  are  satisfied.   In  other  words,  there 


ko  - 


exists  a  function  possessing  the  properties  of  the  fimction  M(e)  associated 
vlth  some  operator  of  type  L  having  m  discrete  eigenvalues. 

Let  us  first  show  the  \iniqueness  of  the  function  M(s).   Suppose  there 
exist  tvo  f-unctions  M  (s)  and  Mp(s),  analytic  and  bounded  in  the  upper  half- 
plane  T  >  0,  having  there  a  finite  number  of  zeros  at  s  =  i/c  and  such  that 

n     n 

M  (-k)   M  (-k) 
^W  =  -^^   =  ^^  ,  (6.10) 

on  the  real  axis;  S(k)  is  a  given  function  satisfying  the  conditions  of  lemma 
6. 5.   Then  the  function  M^(s)/Mp(s)  is  analytic  and  bounded  in  the  upper  half- 
plane  and  real  on  the  real  axis  so  that  it  has  a  bounded  analytic  continuation 
Into  the  lower  halfplane.   By  Liouville's  theorem,  it  follows  that 
M,  (s)/Mp(s)  =  C.   Thus,  the  asymptotic  behavior  of  M,  and  M^  as  |  s  |  -»■  oo 
implies  that  C  =  1,  i.e.  M^  =  Mp.   By  an  analogous  argument,  one  may  prove  the 
uniqueness  of  a  function  M(s)  satisfying  (6.1)  for  given  W(k) .   It  is  only 
necessary  to  consider  ln(M  /Mp). 

We  now  complete  the  solution  of  the  stated  problem.  We  begin  by  recon- 
structing M(k)  from  S(k)  and  consider  first  the  case  when  m  =  0  in  (6.9).   If 
we  normalize  the  phase  T)(k)  =  (l/2)  In  S(k)  so  that  ti(0)  =  0,  then  ti(»)  =  0 
and 

00  00 

Ti(k)  =  -  /  7(t)sln  kt  dt,    r  |7(t)|dt  <  »  (6.11) 

o  o 

by  the  Wiener-Levi  theorem.   The  function 

CO 

M(k)  =   exp   r  r(t)e^^*  dt  (6.12) 

o 

is  the  solution  of  ovoc  problem.      Suppose  now  that  m  ^  0,    the  function 


- 1+1 


n?=l 

possesses  all  the  properties  of  S(k)  but  for  it  m  =  0.   Therefore  the  rela- 
tion 

s(k)  =  %4  (6.11.) 

M(k) 
holds  where  M(s)  has  no  zeros  In  the  upper  halfplane.   The  solution  of  our 
problem  will  then  be  the  fuQction 

M(k)  =  M(k)  \[        ^-^       .  (6.15) 

n=l 

Consider  now  the  construction  of  M(k)  for  a  given  W(k).   By  the  Wiener- 
Levl  theorem^  the  function  p(k)  =  In  W(k)  can  be  represented  by 

00 

p(k)  =  -2  /   7(t)cos  kt  dt.  (6.l6) 

o 

The  function 

00 

M(s)  =  exp  I     7(t)e^^*dt  (6.1?) 

o 

is  analytic  in  the  upper  half-plane^  has  no  zeros  there,  and  possesses  the 
right  asymptotic  behavior  for  leirge  |s|.  The  solution  of  our  problem  will 
be  the  function 

M(k)  =  M(k)    n     mr  ■  ■  ^^-1^) 

n=l 
In  the  following^  the  functions  M(k),  W(k),  and  S(k),  related  to  some  opera- 


k2 


tor  L,   vlll  be   called  the  M-^   W-^    and  S-functions  of  this  operator.      Any  of 
these  functions   characterize  the  spectrum  of  the   corresponding  operator  1,, 
with  M(k)    giving  the  most  absolute   characterization. 

The   following  assertions   are  a  consequence  of  the   above   cited  derivation. 

1)  If  L  has  no  discrete  spectrum,    its  S-functlon  is  uniquely  determined  by 
its  W-functlon  according  to   (6.I6),    (6.I7),    and   (6.1^). 

2)  If  tvo  operators  L     and  L     have  the  same  W-function  ,   L,    has  no  discrete 

o 
spectrum,   but  Lp  has  discrete  eigenvalues  at  X     =   -   ;c^   (n  =  1,  o..,m),    then 

their  M  and  S-functlons  are  related  by 


""   k  +  i/<  ^  /k  -  iK  \  2 


M. 


iW-M^W  n    y-z^,  s^w.s.w  n  hnnr     •  (6.19) 


11=1  n=l  ^ 


These  statements  will  be  used  in  Sees.  8,  11,  and  12. 


7.   The  transformation  operator  A(x,y). 

Together  with  U,,,  in  Sec.  k   there  was  introduced  an  operator  V  with  a 

kernel  in  the  x-representation  given  by  (^.lO).   In  contrast  to  U^,  the 

B 

operator  V  is  not  a  transformation  operator  in  the  sense  of  Sec.  5«   Iii 

B 

fact,  its  kernel  can  not  be  obtained  as  the  Fourier  sine  transform  of  a 

solution  of  (1.1)  vhich  satisfies  a  zero  boundary  condition.   The  present 

section  is  devoted  to  a  clarification  of  the  relationship  betveen  U^^  and  V^. 

B       B 

We  will  suppose  that  there  is  no  discrete  spectrum.   This  allows  one  to  use 

the   concise   operator  notation  vhich  was  introduced  in  the  preceding  sec- 

ti  ons . 

Consider  first  the  operator  U^,  =  U_,W  vhich  is  the  transfonnation  operator 

B     B 

related  to  the  solution 


-  k3  - 


9(x,k)  =  q)(x,k)W(k)   =  2^ 


f(x,k)       f(x,-k) 
M(k)  M(-k) 


(7.1) 


/N^  ''^J 


The  kernel  of  the  operator  K  =  U     -  I  in  the  x-representatlon  is   obtained  by 

B 

taking  the  Foxirler  transform  of  $(x^k)    -  k~     sin  kx: 


K(x,y)  =  I 


J      (?(xA)  -  ^i^)  ^iiL^  k^dk 


art 


f(x,k)    IkxN  -iky^, 


(7.2) 


isx, 


By  virtue  of  the  analyticity  and  boundedness  of  g(x^s)  =  f(x,k)e~   /m(s)  In 

the  upper  halfplane,  K(x,y)  =  0  for  x  >  y.   Therefore,  in  the  x-representa- 

tion  the  operator  U^  is  given  by 

B 


U^f(x)   =  f(x)   +     /     K(x,y)f(y)dy. 


(7.3) 


The   conrpleteness   and  orthogonality  relations 


UgWU*  =   I 


(7.M 


in  terms  of  U  ajre  then  given  by 


UgVro|  =  1,    w  =  w"^ 


(7.5) 


and  from  this. 


(u|)-Vub)"^  =  I. 


(7.6) 


kk 


We  note  that  the  operator  (U*)~  is  a  Volterra  operator  similar  to  U_, 
I.e.  the  integration  goes  from  0  to  x.   In  Sec.  9,    it  will  he  shown^  that 
a  Volterra-type  operator  is  ixniquely  determined  hy  (7.6).   On  the  basis  of 
thls^  one  can  conclude,  by  comparing  (7.^)  and  (7.6)  that 

(U|)~^  =  U^.  (7.7) 

This  establishes  the  relationship  between  the  tvo  operators  U  and  U^.   On 

B      B 

the  other  hand,  by  taking  the  inverse  Fourier  transform  of 

00 

X 

one  easily  gets  the  relation  between  U„  and  V  .  Let  Il(t)  denote  the  function 

B       B 

00 

-(*'  -^  Jim-') ''"'"^-  f^-^' 

—  00 

Because  i/m(s)   Is   analytic  in  the  upper  halfplane, 

n(t)   =0,      t  <  0.  (7-10) 

Recalling  that 

00 

f(xA)    ^  e^^  +    J    A(x,j)e^^^6.j,  (7.II) 

X 

■we  obtain  from  (7.8)  on  the  basis  of  the  convolution  theorem 

y 

A(x,y)  +  n(y-x)  +  J   A(x,t)n(y-t)dt  =  K(x,y).  (7-12) 

X 

This  is  the  sought  relation.   It  will  be  convenient  to  write  this  equation 
as  an  operator  equation.  With  this  in  mind,  we  associate  with  n(t)  the 
operator  Q  which  is  given  in  the  coordinate  representation  by: 


1^5  - 


Qf(x)   =   (I  +  n)f(x)   =  f(x)   +      r    n(y-x)f(y)dy.  (7.15) 

X 

By  virtue  of  the  boimdedness  of  i/m(s),    It  is  not  difficult  to  sIlow  that 
Q=I+II  is  a  bounded  operator  in  our  Hilbert  space.      The  veil-known  properties 
of  integral  equations  -with  difference-type  Volterra  kernels  imply  that  Q  has 
an  inverse 

P  =  I  +  r  =  Q"^  (7.1^) 

vhlch  in  the  x-representatlon  has  a  structure  similar  to  Q: 

00 

Pf(x)   =   (I  +  r)f(x)    Hf(x)   +    Jr{Y  -  x)f(y)dy.  (7. 15) 

X 

The  function  r(t) 

00 

r(t)   -^   f     (M(k)    -   l)e-^^*dk-  (.7.16) 

-  00 

was  introduced  in  the  previous  secticn.   In  terms  of  the  operator  Q  equation 
(7.12)  takes  the  form 

Ug  =  VgQ.  (7.17) 

The  completeness  and  orthogonality  of  the  eigenfunctions  can  now  be  written 
down  using  the  operator  V„.   Substituting  (7-17)  in  (7'5)^  we  find 

VgQW"Vv*  =  1.  (7.18) 

Let  us  clail-fy  the  structure  of  the  operator  QW  Q*.  For  thla,  the  x-repre- 
sentation  will  be  convenient.   In  this  representation,  the  operator  W   =  W 

'">>  jn^  f^ 

has  the  form  W  =  I  +  f]  where  fl  is  an  integral  operator  with  the  kernel 


k6  - 


n(x,y)  .  I  J  2i^  (^^  .  l)  £i^  ,2^  „  H(x-y)  -  H(xH.y)       (7.19) 


and 


-  00 

Denote  by  H^^  the  operator  vlth  kernel  H(x-y)  and  by  H  the  operator  with 
kernel  -H(x+y) .   Then^,  In  our  terms^  the  identity 

w(syMrky  =  ^(-^^  ('^•^i) 

can  be  vritten  in  the  form: 

(I  +  n)(l  +  H^)  =  I  +  r*.  (7.22) 

From  this_, 


QW"-^Q*  =  (I  +  n)(l  +  Hj_  +  H2)(l  +  Ti*)   = 

=  (I  +  r*  +  (i+  n)H2)(i  +  n*)  =  I  +  (I  +  n)H2(i  +  n*).      (7-25) 

One  can  easily  verify  that  the  second  term  in  (7.25)  is  an  integral  operator 
vhose  kernel  is  a  sum  which  can  be  constructed  using  the  function 


F(t)=i  r   (s(k)  --^)e^*dk.  (7.2i+) 

^"^    J        ^  M^(k)  ' 

For  t  >  0^  F(t)  coincides  with  the  function  F(t)  introduced  in  the  preceding 
section.  We  will  associate  with  the  latter  the  operator 


Ff(x)  =  r  F(x  +  y)f(y)dy.  (7.25) 


ifT  - 


We  have  thus  shown  that  completeness  and  orthogonality  of  the  elgenfunctions 

can  he  written  in  terms  of  V  in  the  following  way: 

a 

V^(I  -  F)V^  =  I.  (7.26) 

To  conclude  this  section,  we  note  the  significance  of  the  various  opera- 
tors used  in  Hllbert  space.   By  virtue  of  the  uniform  boundedness  of  M(k), 
l/M(k),  W(k),  l/w(k),  and  S(k),  the  operators  M,  W,  S,  W~  ,  m"  ,  Q,  P,  and 
I-F  are  hounded  in  our  Hllhert  space.  All  the  remaining  operators  used^  to 

wit  U„,  U^,  and  V^,  can  be  obtained  by  multiplying  the  unitary  operator  U^ 
B   B       B 

or  U^''"'^  by  one  of  the  above-mentioned  operators.   Conseguently,  all  of  these 
operators  axe   also  bounded  operators  in  our  Hllbert  space. 


8.   Integral  equations  for  the  kernels  K(x,y)  and  A(x,y). 

In  this  section,  we  shall  establish  how  the  functions  W(k)  and  S(k),  as 

well  as  the  kernels  K(x,y)  and  A(x,y),  are  connected.   The  completeness  and 

orthogonality  relations  for  the  elgenfunctions  written  in  terms  of  U  and 

V  are  such  expressions.  However,  if  we  write  out  th,eBe  equations 
B 

in  the  x-representation,  for  example,  the  kernels  K(x,y)  or  A(x,y)  will 
enter  non-linearly.   This  turns  out  to  be  Inconvenient  in  solving  the  in- 
verse problem.  Nevertheless,  the  fact  that  the  operators  U  and  V„  are 
Volterra  operators  in  the  x-representation  allows  one  to  easily  obtain  ex- 
pressions relating  W(k)  and  S(k)  and  the  kernels  K(x,y)  and  A(x^y)  in  which 
the  latter  enter  linearly. 

For  simplicity,  it  is  first  assxaned  that  there  is  no  discrete  spectrum. 
Rewrite  the  equality 

U^WU^  =  I  (8.1) 

n      n 


-  kQ 


in  the  form 

By  expressing  (8.2)  In  the  x- representation,  one  obtains 

X 

K(x,y)  +  n(x,y)  +   /  K(x,  t)n(t,y)dt  =^K(x,y).  (8.3) 

o 

00 

n(x.y)  =  f  /  ^^  (W(k)  -  1)  ^i£^  k^dk.  (8.M 

o 
Since  K(x,y)  =  0  for  x  >  y,  it  finally  follovs  that 


.A. 

K(x,y)  +  a(x,y)  +   /  K(x,t)n(t,y)dt  =0,   x  >  y.  (8.5) 


In  an  analogous  vay_,  the  relation 

Vg(l  -  F)V*  =  I  (8.6) 

gives 

00 
A(x,y)  =  F(x+y)  +  J   A(x,t)F(t+y)dt,  x  <  y,  (8.7) 

X 

00 
F(t)  =  ^  /  (S(k)  -  l)e^^*dk.  (8.8) 

—  00 

In  case  there  is  a  discrete  spectrum,  similar  equations  can  be  gotten 
by  starting  directly  from  the  representation  (^.5)  ^or   cp(x,k)  and  the  com- 
pleteness and  orthogonality  relations  in  the  form  (2.6).  Equations  for  the 
transfoimatlon  operator  which  relates  the  solutions  cp(x,k)  of  two  different 
equations  of  type  (l.l)  will  now  be  derived.   The  initial  equations  are  the 


k9 


relations { 


.A. 

cpgCx.k)  =  cpi(x,k)   +     /     K(x,t)cp^(t,k)dt, 


(8.9) 


X 


cp^(y,k)   =  cPgCyA)   +     /     K(t,y)cp2(t,k)dt, 


(8.10) 


m. 


XI     C^  ^n  ^''^^n  ^^^  "^  f     f    cp^(x,k)W^(k)cp^(y,k)k^dJi  =  6(x-y 
n^=l       i      i  1  ^ 


)      (1  =  1.2), 


(8.11) 


Multiplication  of  (8.10)  by  cp   (x,k)W„(k)    and  integration  vith  respect  to 


'2^^'^'"2 


k  gives 


m.^ 


2 

■Si 


cp2(x,k)W2(k)cp-L(y.k)k^dk  +     2-j       ^q  ^^(x.iK^   )cp^   (y)   = 


n2=l         2  •"  "2     "2 


A. 

5(x-y)   +      /     K(t,y)6(x-t)dt  =  5(x-y),      x  >  y. 


(8.12) 


Analogoiis   operations  vlth  (8.9)   nov  yield 


m^ 


5(^-y)  =  -     /     cp  (x,k)W  (k)cp  (y,k)k^dk  +    ^     C     cp  (x,1k     )cp  (y,iK     )   + 
^  n^=l     ^2  -^  ^2     ^  ^2 


+    j     K(x,t) 
o 


I      f    9i(t,k)W2(k)cp^(y,k)k^dk  + 


+    2    c^  9l(t,l/c^  )cPi(y.i/c^  ) 


n^^l       2 


dt^      X  >  y. 


(8.13) 


50  - 


Subtracting   (8.11)    for  1=1  from  (8.15)^   ve  obtain 

00 

K(x,y)   +  n(x,y)   +      /     K(x,t)n(t,y)dt  =0,      x  >  y,  (8.1^) 

o 

where 


00 


2 


2. 


i^(x,y)    =  -    j      ^^(xA)    W^Ck)    -  W^(k)   cp^(y,k)k''dk  + 
o 

ng^l      "^2    1  ^2      ^  ^2  ^1      ''l   ''l  ''l 

Since  the  generalization  of  (8. 7)  vlll  not  be  reqizlred  later  on  Its 
derivation  is  omitted.   Equations  (8.5)  and  (8.1^)  are  called  the  Gelfand- 
Levitan  equations,  and  (8.7)  is  called  the  Marchenko  equation. 

These  equations  enable  one  to  solve  Inverse  problems,  i.e.  reconstruct 
the  operator  L  in  different  fonmilatlons.  We  first  apply  these  equations  in 
solving  the  following  two  problems. 

1)  Given  th.e  function  W(k)  possessing  the  properties  enumerated  in  lemma  6.2. 
To  construct  an  operator  L  not  having  a  discrete  spectrum  for  which  W(k)  is 
its  W-fuQctlon. 

2)  Given  an  operator  L,  with  no  discrete  spectrum,  its  W-functlon  and  the 

function  q-,(x)  ^0.  To  construct  an  operator  L  whose  W- function  coincides  with 

2 
that  of  L^  and  which  has  discrete  eigenvalues  at  m  given  points  X     =  -k 

(n  =  1, ...,m). 

The  basic  problem,  i.e.  reconstruct  L  from  its  S-functlon  and  its  dis- 
crete energy  levels,  will  be  solvable  if  these  two  problems  can  be  solved. 
For,  let  there  be  given  a  function  S(k)  possessing  the  properties  enumerated 
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in  lenma  6.3  and  m  distinct  positive  nimbers  k     (n  =  !_,..  .^m) .      By  the  pro- 
cedure described  in  Sec.    6,    a  unique  W(k.)    can  be   constructed  from  this  data. 
If  in  conjuziction  with  problem  1^   we  now  construct   the  operator  L  without  a 
discrete  spectrum  for  which  this  W(k)   is  its  W-functlon  then  by  the  asser- 
tions of  Sec.    6,    the  S-function  of  this   operator  will  be 

^  /k  -  iK  \2 

^^^)  =  ^^^^  n(kTi^  •  (8.16) 

If  in  accordance  with  problem  2,   we  now  start  with  the  operator  L-,  and  construct 

the  operator  L  which  has  the  same  W-function  as  L,  and  m  discrete  eigenvalues 

2 
at  the  points  \     =  -k  ,    then  L  will  have  as  its  S-function  the  initially-given 

function  S(k)  and  the  basic  problem  will,  therefore,  be  solved. 


^  9 •   Existence  of  a  solution  of  the  Gelfand-Levltan  equation.   Solution  of  the 
first  problem. 

Let  us  proceed  to  solve  the  problems  formulated  in  the  previous  section. 
To  do  this,  we  will  use  the  Gelfand-Levltan  equations  (8.5),  (8.4)  and  (8.1^), 
(8.15)  and  we  therefore  first  ahow  the  existence  of  solutions. 

We  go  directly  to  the  general  equation- (8. l4) .   Let  the  operator  L^  be 

given.   That  is,  its  eigenfunctions  cp, (x,k),  its  W-functlon,  W  (k).  Its  dis- 

2 
Crete  eigenvalues  X  =  -K       and  the  corresponding  normalizing  factors  C 

^1     ^1  ^1 

are  given.  Also,  let  the  function  W(k)  possessing  the  properties  enumerated 

in  lemma  6.2  as  well  as  the  arbitrary  positive  numbers  k     and  C  be  given. 

Assume  that  the  k     are  distinct.   On  the  basis  of  this  data,  we  construct  the 

n  ' 

function  f2(x,y)  by  means  of  (8.I5).   By  virtue  of  the  given  conditions,  this 
function  will  be  absolutely  Integrable  with  respect  to  x  or  y  in  any  finite 
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Interval.   Equation  (8.1^)  for  K(x^y)  Is  an  equation  with,  respect  to  the 
argument  j,    the  kernel  and  free  term  depending  on  x  as  a  parameter.   For 
fixed  x^  this  equation  is  a  Fredholm  equation.   Therefore^  to  shov  that  a 
unique  solution  exists,,  it  suffices  to  prove  that  the  homogeneous  equation 
has  only  a  trivial  solution. 

Suppose  that  for  fixed  x  ,  the  equation 


h^(y)  +  /   h^(t)fl(t,y)dt  =  0 


(9.1) 


has  a  solution.  Denote 

h^(y)   for  y  S  x^ 

ii(y)  = 

for  y  >  X 


(9.2) 


Multiplying   (9.1)  by  h(y)   and  integrating  with  respect  to  j,   ve  obtain 


r    li^(y)(iy  +     jj    li(t)n(t,y)h(y)dtdy  =  0, 


o  o 


(9.5) 


or,  if  ve  substitute  the  expression  for  fi(x,y)  from  (8. 15) 


J    h^(y)dy  +  ~  J        J    h(y)cp^(y,k)dy 


•-  o 


W2(k)  -  W^(k) 


k^dk  + 


^ 


m      p  0°  n  p     -L       00 

+  2  ^n  /  li(y)^i(y^l'<n^'^^   ■  -2  CJn  J   ^^^^^n  ^^^"^^ 


=  0. 


(9A) 


By  the  application  of  Parseval's  equality  for  the  system  of  fimctlonB  cp^(x,k) 
ttde  may  be  written  in  the  form 
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OOpOO  — r  mpoo  - 

oLo  J  °^lLo 


=  0. 
(9.5) 


Since  Wp(k)  Is  positive^  this  Implies  that  h(y)  is  orthogonal  to  the  sub- 
space  corresponding  to  the  continuous  spectrum  of  L.   Consequently,  Wp(k) 
is  a  linear  combination  of  eigenfunctions  of  the  discrete  spectrum.  But  this 
is  impossible  since  h(y)  vanishes  identically  for  j  >  x  ,    so  that  h(y)  =  0. 
The  rem.ainder  of  this  section  and  the  tvo  succeeding  sections  will  be 
devoted  to  solving  the  first  problem..   In  this,  we  will  be  dealing  with  the 
equation 


K(x,y)  +  fi(x,y)  +  /  K(x,t)r2(t,y)dt  =  0,  x  >  y. 


{9.6) 


^/    \   2   r  sin  kx  /,,,,  ^   ^  v  sin  ky  ,  2 
f^(x,y)  =  ^  /  ^-  (W(l^)  -  1)  ^  ^  ^, 


(9.1) 


which  was  obtained  from  the  condition 


u^m*  =  I 


(9.8) 


for  the  operator  U„  =  I  +  K.   Conversely^  we  shall  show  that  the  operator 

n  =  I  +  K  constructed  with  the  help  of  the  solution  of  this  equation, 
B 

possesses  the  property  [9.Q),   where  the  function  W(k)  is  only  required  to 
satisfy  the  conditions  enumerated  in  lemma  6.2. 

We  note  that  the  solution  of  (9'6)  determines  a  kernel  K(x,y)  which 
is  different  from  zero  only  for  x  <  y: 


K(x,y)   =      < 
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X 

K(x,y)   +  n(x,y)   +      /     K(x,t)n(t,y)dt,      x  <  y, 

o 


{9.9) 
0,  X  >y. 


With  tMs  kernel^   we  will  associate  the  operator  K.     Let  U     =  I  +  K.      In 

B 

terms  of  U^,  and  U^,^  relation  (9- 9)  becomes 

U3W  =  U^,  (9-10) 

and  to  prove  (9.8),  it  remains  to  be  shown  that 

Ug  =  (u*)"^  (9.11) 

Since  W(k)  Is  real^  W  is  self adjoint.  From  (9-10)^  it  followB  that 

U„U*  =  U„WU*  so  that  U^U*  is  also  self ad joint.   On  the  other  hand,  since  U^ 
B  B     B   13  B  B  B 

and  U*  are  Yolterra  operators,  U^U*  is  also  of  Yolterra-type  and  this  pro- 

perty  together  with  the  selfadjolntness  Implies  that  U„U*  =  I.   This  is  what 

B  B 

we  had  to  show. 

These  arguments  are  of  a  rather  formal  character  since  no  estimates  for 
the  kernels  K(x,y)  and  K(x, y)  were  given  which  would  have  explained  the  mean- 
ing of  U^  and  U^  as  operators  in  Hilbert  space.  However,  relation  (9«8)  when 

B      B 

written  out  in  the  x-representation  only  Involves  integrals  over  finite 
intervals  of  integration  (whose  convergence  it  is  not  necessary  to  prove). 
And  this  relation  can  be  pi^Dven  starting  from  {9'^)   directly.  Anyway,  inas- 
much as  this  relation  is  proven,  it  is  possible  to  attach  meaning  to  the 
operator  U^.  In  fact,  the  function  M(k)  constructed  from  W(k)  by  the  proce- 
dure  of  Sec.  6,  Is  a  bounded  function  and  a  bounded  operator  can  be  associ- 
ated with  it,  namely,  an  operator  which  multiplies  by  M(k)  in  the  momentum 
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representation.      Since  M(k)   ^  Q,    the  operator  M       is  also  bounded.      From 

''+)  -1 

(9-8)^   ve   find  that  U"'        =  UJ^I       is  unitary  and  from  this   follovB  the  bounded- 

ness  of  U   . 

B 

With  the  operator  thus  obtained^   ve  define  an  operator  L  by 

2 
vhere  L     is  the   operator  vhich  multiplies  by  k     in  the  k-representation  vhich 

vas  introduced  in  Sec.    5  and  Sec.   5.      Since  L  =  U^L  WU*  by  virtue  of  (9-8)^ 

Bo      B 

L  is   selfadjolnt  because  W  and  L     commute.      On  the  basis  of  theorem  5- 2,   ve 
can  assert  that  L  is  a  differential  operator  in  the  x-representation  27elated 
to  an  equation  of  type   (l.l)   vith  a  potential  given  by 

q(x)   =   2^K(x,x).  (9.15) 

Equation  (9.8)  shovs  that  the  given  function  W(k)  is  the  W- function  for  L. 
And  so^  the  first  problem  of  Sec.  8  is  solved  in  the  sense  that  an  operator 
L  has  been  shovn  to  exist.   The  properties  of  the  potential  q(x)  occurring 
in  the  solution  of  this  problem  vill  be  investigated  in  the  folloving  section. 

In  conclusion^  ve  note  that  our  considerations  imply  that  (9-8)  deter- 
mines a.  unique  Volterra  operator  U_,.  This  fact  vas  previously  used  in  Sec.  7« 

B 


10.  Marchenko's  equation.  The  properties  of  the  potential. 

In  studying  the  properties  of  q(x)^  it  is  convenient  to  use  Marchenko's 
equation 

A(x,y)  =  F(x+y)  +  [    A(x,t)F(t+y)dt,  x  <  y,  (lO.l) 
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F(^)  =h     f    (S(l^)  -  De^^V.  (10.2) 


2n 


The  existence  of  a  solution  of  (lO.l)  will  not  be  proven.   Instead^  it  -will 
be  shown  that  It  Is  equivalent  to  the  Gelfand-Levltan  equation  (8.5)  on 
the  condition^  of  course^  that  the  corresponding  functions  W(k)  and  S(k)  are 
related  by  the  formulas  of  Sec.  6.   In  the  preceding  section^  the  equiva- 
lence of  the  Gelfand-Levltan  equation  (8.5)  and  the  relation  U^WU*  =  1  for 

B  B 

the  operator  U^  =  I  +  K  was  shown.   By  repeating  the  arguments  of  Sec.  7 

~  -1 
and  Sec.  8,  one  will  find  that  V  =  U^Q   Is  an  operator  which  can  be  con- 

B  B 

structed  using  the  solution  A(x,y)  of  (10. 1).   Therefore  (10. 1)  has  a  solution 

for  X  g  0.   Conversely,  it  Is  easily  shown  that  if  the  operator  V^  is  formed 

B 

with  respect  to  a  solution  of  (10. 1),  then  U^  =  (Q*V*)~  is  the  operator  form- 

B  B 

ed  by  using  the  solution  of  the  CJelfand-Levltan  equation  (8.5).   This  implies 
the  uniqueness  of  a  solution  of  (10. 1). 

The  potential  q(x)  =  2  —  K(x,x)  can  be  simply  expressed  in  terms  of 
A(x,y).   Thus,  from  (7-12)  It  follws  that  A(x,x)  =  K(x,x)  -  n(o)  and 
dA(x,x)/dx  =dK(x,x)/dx.   If  (7-7)  is  written  out  in  terms  of  K(x,y)  and  K(x,y), 
one  obtains 

K(y,x)  =  K(x,y)  +  J   K(t,x)K(t,y)dt  =  0,  (IO.5) 

y 

from  which  it  follows  that  dK:(x,x)/dx  =  -dK(x,x)/dx.  Finally, 

q(x)  =  2  ^K(x,x)  =  -2^A(x,x).  (IO.I+) 

The  solution  of  (lO.l)  will  now  be  investigated.  It  turns  out  that  q(x) 
behaves  in  many  respects  like  the  derivative  of  F(t)  for  t  >  0.  We  obtain 
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some  estimates  and  first  consider  how  F'(t)  depends  on  q(x).   It  Is  conveni- 
ent to  Introduce  the  fimctlons 


W  BO 

cj(x)  =   /   |(i(t)|dt,   a^(x)  =   /   t|q(t)|dt. 


(10.5) 


The  Inequalities   cited  in  Section  k   for  A(x,y)  and  Its  partial  can  now  be 
written: 


lA(x,y)|  gKa(^) 


(10.6) 
(10.7) 


In  (lO.l);,  set  X  =  y.   Then 


or 


A(x,x)  =-F(2x)  +   /  A(x^t)F(t+x)dt 

X 

00 

F(2x)  =  A(x,x)  -  2  /  A(x^2t-x)F(2t)dt. 


(10.8) 


(10.9) 


The  last  equation  is  solved  by  the  method  of  successive  approximations^  and 
an  estimate  for  its  solution  can  be  gotten 

|f(2x)1  g  Kff(x).  (10.10) 

Differentiation  of  (10.9)  with  respect  to  x  and  the  inequalities  (lO.T)  and 
( 10.10)  then  imply 


F'(2x)  +  J  q(x) 


Ka  (x) 


(10.11) 
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We  now  study  the  converse  question  for  which  It  Is  convenient  to  transcribe 
(lO.l)  to  the  form: 

CO 

A(x,x+y)  =  F(2x+y)  +   /  A(x,x+t)F(2x+t+y)dt.  (10.12) 

o 

Let  F  denote  the  operator 

00 

F^g(y)  =   /  g(t)F(2x+t+y)dt.  (10.15) 


Because   (lO.l)   has  a  iinlque  solution  for  all  x  and  since  F     is  small  for 
large  x^    the  noim  of  (l   -  F   )        In    ^^(0^«>)   is  uniformly  boimded.      This  im- 
plies: 

00  00  00 

J    |A(x,x+y)|dy  =    J     |A(x,y)|dy  ^  K    J       |F(t)|dt.  (I0.l4) 

o  X  2x 

Introducing  now  the  notation 

00  oo 

t(x)  =  J     |F'(t)|dt,   T^(x)  =  J    t|F'(t)ldt,  (10.15) 

X  X 

one  may  write 

CO 

|f(x)|  g  I      lF'(t)|dt  =  t(x).  (10.16) 

X 

From  (lO.l)  and  the  inequality  (lO.l^),  the  ixniform  estimate 

|A(x,y)|  S  KT(x+y)  (lO.l?) 

is  gotten.  Let  us  consider  the  differentiability  of  A(x,x).   Introduce  the 
notation: 
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B(x,y)  =  A(x,x+y) 


(10.18) 


The  difference  quotient 


A^B(x,y)   B(x+h,y)  -  B(x,y) 


h 


(10.19) 


satisfies  an  equation  of  the  same  kind  as  B(x,y) 


A  ^B(x,y)   A  ^F(2x+y) 


+   /   B(x,t) 
o 


A  F(2x+t+y) 


dt  + 


A  B(x,t) 


X 


F(2x+2h+t+y)dt. 


(10.20) 


The  free  term 


A  P(2x+y) 


B(x,y) 


A  F(2x+t+y) 


dt 


(10.21) 


can  he  estimated,  uniformly  in  h,  hecause  of  the  differentiability  of  F(t), 

vhence  follows  the  differentiability  of  B(x,y).  Moreover,  in  exactly  the 

dB 
same  way  as  before,  one  obtains  the  estimate  for  ^(x,y): 


lB;(x,y)|dy  g  Kt(2x), 


(10.22) 


This  implies  that 


^A(x,x)  -  2F'(2x) 


g  Kt  (2x) 


(10.25) 


or 


P'(2x)  +  J  g(x) 


g  Kt  (2x). 


(10.24) 
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Thus   the  inequalities  (10. ll)  and  (10.2^)  explicitly  show  the  analogous 
behavior  of  q.(x)  and  F'(2x).   In  particular,  they  imply  that  if 

00 

r    x|q(x)ldx  <  00  (10.25) 

o 

holds,  then  so  does 

00 

r  x|F'(2x)|dx  <  CO  (10.26), 

o 

and  conversely.   These  inequalities  also  can  he  used  to  establish  necessary 
and  sufficient  conditions  on  the  S-function  guaranteeing  that  the  corre- 
sponding potential  (Je cay  like  x~  ,  or  exponentially,  or  vanish  identically 
for  X  >  A,  etc.  A  precise  formulation  of  these  will  not  be  given  here 
because  of  the  lack  of  space. 

11.  Krein's  equation.  Asymptotic  behavior  of  cp(x,k) 

The  eigenfunctions  cp(x,k)  of  the  operator  L  are  formed  using  the  formiaa 

X 

,  ,  .   sin  kx    P    -u-/        \    sin  ky  m  i^ 

cp(x,k)  =  ^—  +   /  K(x,y)  ^—  dy.  (11-lJ 

o 
By  (1.21),    these  functions  behave  asymptotically  like 

cp(x,k)  «  ^  sln(kx  -  Ti(k)),  (11.2) 

for  large  x.  A(k)  and  T](k)  are  respectively  the  modulus  and  argument  of  the 
M-functlon  of  L  which  is  luaiquely  determined  by  W(k)  by  means  of  the  relation 


61 


and  the  condition  of  analyticity  (see  Sec.  6). 

It  is  not  easy  to  get  the  asymptotic  expression  (ll.2)  directly  from 
(11. l).  Another  representation  of  9(x^k)  will  therefore  he  obtained  from 
which  (11.2)  can  be  deduced  without  difficulty.   Of  course,  this  represen- 
tation will  not  only  serve   to  prove  (ll.2),  but  the  procedure  outlined  be- 
low is  independently  interesting  in  that  it  reveals  more  completely  the 
structure  of  the  kernel  A(x,y)  of  the  transformation  operator. 

The  starting  point  is  the  Gelfand-Levitan  equation 

X 

K(x,y)  +  n(x,y)  +  /  K(x,t)n(t,y)dt  =  0,  y  <  x,  (H.^) 


fi(x,y)  =  -  /  sin  kx[w(k)  -  l]  sin  ky  dk.  (11-5) 


o 


From  (11.5),  it  is  evident  that 

fl(x,y)  =  H(x-y)  -  H(x+y),  (11. 6) 

where 

00 

jj(t)  =  -     f    [w(k)  -  ijcos  kt  dk  (11-7) 

o 

coincides  with  the  function  E(t)  introduced  in  Sec.  6  (see  (6.6)).   If  a 
solution  K(x,y)  is  sought  in  the  form 

K(x,y)  =  r^{x-j)    -  T^{x+y),  (11.8) 

then  one  obtains  for  T^{t)   the  equation 

2x 

r2^(t)  +  H(t)  +  J    r2^(s)H(s-t)ds  =  0,  (11.9) 
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which  Is  equivalent  to  (11.4).   This  Is  an  equation  for  ?_  (t)  In  the  argu- 
ment t,    X  playing  the  role  of  a  parameter.   This  fact  underscores  the  unsym- 
metric  dependence  of  F  (t)  on  its  arguments.   Equation  (II.9)  is  caJJLed 
Krein's  equation.  Knowing  its  solution^  one  can  easily  find  the  kernel 
K(x,y)^  and  together  with  it,  the  potential  q(x)  and  the  solution  cp(x,k). 
However,  these  functions  can  be  expressed  directly  in  terms  of  Pp  (t),  as 
follows . 

If  the  resolvent  of  the  kernel  H(x-t)  on  the  Interval  (0,a)  is  denoted 
hy  G  (x,t): 


a 


then 


CL 

G  (x,t)  +  H(x-t)  +   /   G  (x,s)H(s-t)ds  =  0,  (ll.lO) 

a  J    a 


r^{t)   =   G2^(0,t).  (11.11) 


We  now  note  some  properties  of  G  (x,t).   Differentiation  of  (ll.lO)  with 
respect  to  a  shows  that  G  (x, t)  satisfies 

4-  G  (x,t)  =  G  (x,a)G  (a,t).  (11.12) 

Also, 

G  (t,s)  =  G  (a-t,a-s).  (11.15) 

In  terms  of  rp^('t)  =  G„  (0,t),  q(x)  =  2  dK(x,x)/dx  is  given  by 

q(x)  =  2  ^  G2^(0,0)  -  2  ^  G2^(0,2x).  (ll.l4) 

If  the  notation 

A(x)  =  2G2^(0,2x)  =  2r2^(2x)  (11.15) 
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is   introduced,    then  by  virtue  of  (11.12)    and   (11.15) 


<l(x)   =   -  ~A{x)   +  a2(x) 


(11.16) 


This  last  formula  can  he  used  to  reduce  the  second  order  differential  equa- 
tion containing  q(x)  to  a  system  of  equations  of  the  first  order.   Thus^ 


dx 


-.(x)=  (i-A(x))  (A.A(X)) 


(11.17) 


and  the  equation 


-y"  +  q(x)y  =  k  y 


(11.18) 


is  equivalent  to  the  system: 


^  y  +  Ay  =  kz, 


> 


dx 


z  +  Az  =  ky. 


(11.19) 


In  certain   cases,    this   system  turns   out  to  be  more   convenient  than  the   origin- 
al equation   (II.I8).      Let  us  now  express  the  solution  cp(x,k)   in  terms   of 


x 


/      ,  V        sin  kx  /     „     /        N    sin  ky 

cp(x,k)  =  — ^^  +     /     r2^(x-y)  —^ 


^y  -  /  r^i-^j)  ^^  ciy  = 


2x 

£l£_kx  +      r      r     (t)   "^^  M,x-t)   ^t  = 
k  J  2x^  k 

o 


»i- 


2x 
e'^^fl-.     [      r2^(t)e-i^\t 


(11.20) 


6k  - 


As  X  -*  00,  ve  find 


cp(x,k) 


ii. 


Ikx  /, 
e     1  + 


r(t)e    dt 


(11.21) 


vhere  Vit)   =   lim  T^   (t)   is  the  solution  of  the  equation 

X  -►  00 

00 

r(t)  +  H(t)  +  /  r(s)E(t-s)ds  =  0. 


(11.22) 


That  Is^    this   equation  is   satisfied  by  the   function 


r(t)  =  ^    /     (M(k) 


-   l)e  dk, 


(11.23) 


introduced  in  Sec.  6  (see  (6.5)).   This  is  not  difficult  to  show  if  one  takes 
the  Fourier  transform  of  the  identity  (11.5)  written  out  in  the  form 


(W(k)-l)(M(k)-l)  +  (M(k)-l)  +  (W(k)-l)  =  j;j^  -  1, 


(ll.2i^) 


and  takes  into  account  that  l/M(-k)  is  analytic  in  the  lower  halfplane.   Thus, 
from  (11.21),  it  follows  that 


q,(x,k)  -  i  im  |M(-k)e-^^] 
as  X  -♦  00  and  this,  in  turn.  Implies  the  asymptotic  behavior  (11.2). 


(11.25) 


12.   Relationship  of  the  operators  to  the  discrete  spectra. 

In  line  with  our  program,  we  have  completed  the  solution  of  the  first  prob- 
lem and  in  this  section  we  begin  to  consider  the  second  problem.   The  Gelfand- 
Levitan  equation  will  be  considered  for  the  case  when  the  given  operator  L,  does 
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not  have  a  discrete  spectrum^  i.e.  all  C  =  0,  and  the  operator  being  sought 

has  the  same  function  W(k)  as  L^  and  discrete  eigenvalues  at  X     =   -k 

1  n     n 

(n  =  l^...,m)  with  corresponding  normalizing  factors  C  .  The  equation  in  this 
case  Is  given  by: 

X 

K(x,y)  +  fi(x,y)  +  /  K(x,t)n(t,y)dt  =0,  x  >  y,  (12.1) 

o 

m 

fi(x,y)  =  ^  C^cp^(x,i/<^)cp^(y,iK^).  (12.2) 

n=l 


This  is  an  equation  vlth  a  degenerate  kernel  and  is  easily  solved.  For  thls^ 
it  is  advantageous  to  use  vector  notation.   Thus  we  write  n(x,y)  in  the  form: 

n(x.y)  =  (^(x),9(y))  (12.5) 

where  ijf(x)  and  cp(y)  are  vectors  with  the  components  C  cp  (x^Ik:  )  and  cp  (y^i/c  ) 
respectively.  We  will  seek  K(x^y)  in  the  form 

K(x,y)  =  (a(x),^(y)).  (12.1^) 

Then  the  Gelfand-Levitan  equation  transforms  to 

(a(x),$(y))+(7(x),^(y))+(   y    R(t)dta(x),$(y)  )     =0.  (12.5) 

Here  R(t)  is  the  tensor  product  of  the  vectors  cp(t)  and  ^{t) ,    i.e.  the  matrix 
with  elements 

r^jj.(t)  =  cp^(t)i|r^(t)   (i,k  =  l,...,m).  (12.6) 


6G 


Let 

X 

V(x)   =  I  +      /     R(t)dt.  (12.7) 

o 

[phen  the   general  result   concerning  the  existence   of  a  solution  of  the   Gel- 
fand-Levitan  equation  implies  that  the  matrix  V(x)    has   an  inverse   for  all 
X.      Moreover,    from  its   definition,    one    can  easily  verify  that  this  matrix 
is  positive-definite. 

By  virtue   of   {V2..k)    and   (12.5),   ve   find 

K(x,y)   =   -(V(x)"^?(x),^(y)),  (12.8) 

and  in  particular. 


K 


(x,x)   =   -(V(x)~    ^(x),9(x))   =   -   sp(v(x)~^(x))   = 


sp(v(x)'^  ^  V(x))    =   -     ^lndetV(x),  (12.9) 


so  that 


R(x)   =  ^  V(x).  '  (12.10) 

On  the  basis   of  an  Investigation  of  the   general  Gelfand-Levitan  equation 
(8.1^)    analogous   to  that   carried  out  for  (8.5) >    one    could  shov  that  the 
operator 

L  =  UgL^U^^  (12.11) 

formed  by  using  the   solution  K(x,y)    of  this   equation  is   a  differential  opera- 
tor in  the  x-representation  with  a  potential  q(x)   =  q,  (x)   +   Aq(x)   where 

Aq(x)   =   2  ^  K(x,x)       .  (12.12) 
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Moreover,,  its  eigenfunctions  are  given  by 


A, 


(12.15) 


its  W- function  Is  the  function  w(k)  and  the  quantities  X     =  -K^   (n  =  1, . .  ..,m) 
are  the  points  of  its  discrete  spectrum.   Such  an  investigation  will  not  be 
carried  out^  and  the  proof  of  the  above  assertions  for  the  problem  under  con- 
sideration will  be  effected  by  other  means  which  are  more  simple.   The  fact 
that  the  function  cp(x^k)  given  by  (12.13)  is  a  solution  (l.l)  with  a  potential 
q(x)  =  q-,(x)  +  Aq(x)  where  A  q(x)  is  given  by  (12.12)  will  be  deduced  in  Sec. 
15.   In  this  section,  the  properties  of  these  functions  will  be  examined. 

From  (12.8)  and  (12.9),  one  obtains  the  following  formulas  for  Aq(x) 
and  cp(xjk) : 


A  q(x)   =   -   2  -^  In    II  V(x)  ||,    ||  V(x)   ||   =   det  V(x)  , 


dx 


;i2.iM 


cp(x,k) 


V(x)  ?(x) 

p(x,k)      9-L(x,k) 

II  V(x)   II 


(12.15) 


The  numerator  of  this  last  equation  is  the  determinant  of  the  matrix,  composed 
from  the  matrix  V(x),  bordered  by  the  vectors  f(x)  and  p(x,k)  where 


p(x,k)  =  /  ?(y)(P;j_(yA)(iy- 


(12.16) 


The  M- function  of  L  can  be  found  from  the  asynrptotic  expansion  of  the  solution 


68 


9(x,k).   We  will  next  investigate  this  asymptotic  behavior  and  also  study 
the  properties  of  the  potential  Increment  ^q(x)  directly  from  (l2.l4)  and 
(12.15)  restricting  ourselves  for  simplicity  to  the  case  m  =  1.  We  will 
make  use  of  the  fact  that  cp  (x^  la),  a  >  0,  behaves  asymptotically  like 


q)^(x,la)   =  Ne      [l  +  o(l)] 


(12.17) 


as  X  ->  00  ,    and  as  x  -*  0   ^ 


cp^(x,ia)   =  x[l  +  0(1)]. 


(12.18) 


Express  cp(x_,k)  a.s  the  sum  of  two  terms  : 


cp(x,k)  =  2J2  [li(x,k)  -  li(x,-k)]. 


(12.19) 


where 


Ikx 


h(x,k)  =  M(-k)e"^  +  0(1) 


(12.20) 


as  X  -*  »   .      Equation  (12.I5)   then  implies   that 


h(x,k)  =  h^(x,k)   - 


/  cp^(x,lK)h^(x,k)dx  = 


1  +  C   /    [cp^(t,i/<)]^dt 


=  M^(-k) 


ikx 


CNe' 


KX 


N 


lkx+/oc 


„2     ^         <  +  Ik 
^  N       2kx 

2/c 


+  0(1) 


„  ,   ,  N   k  +  l/c     Ikx   .      /^  ^ 


(12.21) 
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BO  that 

M(k)   =  M^(k)   ^^   .  (12.22) 

2 

From  this,    it   follows   that  L  has   a  discrete  eigenvalue   at   the  point   \  =   -<   . 

Let  us  now  consider    Aq(x).      By  virtue  of  its   definition   (12.7),   V(x) 
is  twice  different! ah le  even  if  the  potential  q  (x)   of  the  given  operator  is 
a  generalized  function  with  6-type   singularities.      Therefore    Aq(x)    is   always 
locally  summable.     As  x  -►  0     ,   we  have 


<1  +  C  /    t^[l  +  o(l)]dt>     =   -4Cx[l  +  o(l)]. 
I  o  , 


Aq(x)   =   -2  -^  In     {l  +  C   /   t''[l  +  o(l)]dt^     =   -4Cx[l  +  o(l)].  (12.25) 


and  as  X  -♦  « 


A<,(x)   =   -2  V"(x)V(x)-(V(x))^  ^  _  2   ^^^p  ^-2.x^^  ^  ^^^^^^  ^^^^.^^ 

V^(x)  ^ 

In  the  general  Gase,m  >  1,    asymptotic  forms   similar  to   (12.25)    and  (12.2^)   still 

m 

hold  for  Aq(x)  in  which  C  =  Z.  C  in  (12.25)  and  C  =  C  ,  k:  =  K^   in  (l2.2ij-), 

n=l 

K     heing  the  smallest  of  the  k   .   In  eveiy  case, 
r  n 

00 

r  xlA(i(x)ldx  <  CO  .  (12.25) 


The  solution  of  the  second  prohlem  of  Sec.  8,  except  for  the  one  statement 
to  be  proven  in  Sec.  15,  is  thus  completed.   The  solution  is  given  by  (12.1^) 
and  (12.15).   In  these  formulas,  arbitrary  positive  numbers  may  be  inserted  for 

the  C  so  that  we  have  an  m-parameter  family  of  solutions.  Formula  (I2.l4)  is 

n 

the  formula  for  the  equivalent  potentials  given  by  Jost  and  Kohn. 
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This  Blnrultaneously  conpletes  the  solution  of  the  basic  Inverse  problem. 
The  resultB  obtained  in  Sees.  9-12  are  summarized  in  the  following  theorem. 
Theorem  12.1.  Every  function  S(k)  possessing  the  properties; 

1)  |s(k)l  =  S(»)  =  S(0)  =  1, 

2)  S(-k)  =  S(ky  =  S"^(k), 

00  oo 

5)   S(k)  =  1  +  /  F(t)e"^^'^dt  vlth  /   |F(t)ldt<", 


h)      arg  S(k) 


00   ■•.'» 


=  -  ^irtm,  m  g  0 


is  the  S-function  of  some  operator  having  a  continuous   spectrum  along  the 
half-ray  (O^")    and  m  negative  eigenvalues.      In  the  x-representation^    it  is   a 
differential  operator  of  type    (l.l)   with  a  potential  that  may  be   a  generalized 
function  such  as  the  derivative  of  a  locally  s-ummnble  function. 

00 

In  order  that  the   condition        /     x|q(x)|dx  <  °°  be   fulfilled,    it  is  neces- 

"  o 

sary  and  sufficient  that      /     t|F'(t)|dt<«'.      Ifm>0,    the  potential  is  not 


o 

S£      „  -  .  ....  .    i_. / 

O 

uniquely  detennined.      An  m-parameter  family  of  potentials  exists  such  that  the 


associated  operators  have  S(k)    as  their  S-functlon  and  the  given  quantities 


2 
)t     =  -k:     as  eigenvalues, 
n  n 


15.      Operators  with  M- functions  differing  by  a  rational  factor. 

IBae  general  G^lfand-Levltan     equation  could  be  solved  in  closed  form  for 
the   case  treated  in  the  preceding  section.      This  is  not  the   only  situation  in 
which  this  is  possible.      It  turns  out  that  one  may  explicitly  solve  the  opera- 
tor equation  for  the  transformation  which  connects   two  operators  L^    and  L_ 
whose  M-functions  differ  by  a  rational  factor.     We  will   consider  the   case  when 
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the  given  operator  L^  with  a  potential  q.,(x)  does  not  have  a  discrete  spectrum 
and  we  will  construct  the  operator  Lp  whose  M-functlon  Is  given  hy 


TT   15:  +  la, 


k  +  ip.  ^   '^l>°^  h   >°- 


(15.1) 


.•t=l 


The  requirement  3^  >  0  guarantees  that  M(k)  will  be  regular  In  the  upper 
halfplane.   The  condition  a  >  0  does  not  disturb  the  generality  of  the  dls- 


cusslon  since 


k  -  la  _  k  -  la  k  +  la 
k  +  ip  ~  k  +  la  k  +  13  ^ 


(15.2) 


one  may  consider  the  transformation 


"O')  -  "iC^'  Iflf  ' 


(13.3) 


and  afterwards  carry  out  the  transformation 


•'sO^)  -  "i"^)  mi ' 


(I3.it) 


as  was  done  in  the  preceding  section. 

We  will  also  assume  that  all  the  a,  and  p  are  distinct.   The  equation 

V  V 

for  the  related  transformation  operator  is  given  by: 


K(x,y)  +  li(x,y)  +-.  /  K(x,t)a(t,y)dt  =  0,   y  <  x 


(15.5) 


n(x,y)  =  f  r  q>i(xA)  [W2(k)  -  W^(k)jcpi(y.k)k^dk  = 


=  ^  J    cp^(x,k)Wi(k) 


1=-  1  ^  +  =^1 


<p^(y,k)k  dk. 


(15.6) 
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If  the  term  in  the  brackets  is  resolved  into  partial  fractions: 


vhere 


N       ,  2        „2                     N 

il          2          2                   -<--< 
^^^     k     +a^                    1=1 

?          2  ' 
k     +a^ 

N 

n  of  -  -j) 

n  ^  2      2. 

(13.7) 


(15.8) 


ii(x,y)  can  be  expressed  in  the  form 

n(x,y)  =  f  r  cp^(x,k)  2j   2  ^  2  W^(k)q5j_(y,k)k^dk  =  ]^  A^g^  (x,y).  (15. 
^  t=l  k  +  a„  t=l     t 


9) 


2 
Here^  g  (x,y)  is  the  resolvent  kernel  of  the  operator  L  for  X  =  -a  .   In  Sec. 

2,  it  vas  seen  that 

fj_(y.ia) 
g^(x,y)  =  cp^(x,ia)  -^  ,^^s^     ,  x  <  y,  €^(x,y)  =  g^(y,x)  ,     (15.IO) 


is  a  solution  of 


^  +  a^  +  q^(x)  j    gjx,y)  »  5(x-y).  ^^^'^^ 

Also,  any  solution  of 

-i|r"(x)  -  k%(x)  +  (l^(x)i)r(x)  =  0  (13.12) 
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satisfies 


(a^  +  k^)  J    gjx,y)i|r(y)dy  =  ^(x)  +  [s^(x,y);  t(y)]  |y^^  ,      (15.13) 
a 
vhere 

[cp(x);^(x)]  =  t'(x)lr(x)  -  cp(x)t' (x).  (15.1^+) 

We  will  seek  a  solution  K(x,y)  of  (15-5)  In  the  form: 

N 
K(x,y)   =     2Z  a  (x)cp5(y),  (15-15) 


J=l     ^'      ^J 


vhere 


cpjCy)  =  q'lCy^iPj)    (J  =  i,-...N).  (15.16) 

The  substitution  of  (15-15)   into   (15-5)   gives: 

j=i    -J       J  t=i  mJ^ 

N       N  X 

+     SX]A,a(x)      r    q)P(t)g^(t,y)dt  =  0.  (15.1?) 

J=l   t=l      ''   ^  ^Q        ^ 

By  (15.15),  the  last  term  can  he  transformed  as  follows: 

XIS^.^w  r<p?(t)g?(t,y)<it  =  2-iZ!^(^)  -2-S<p?(y)  + 


J=l  i=l  -^    L         '   -■  ^t  -  Pi 


-  7^ 


The  first  term  in  this  cancels  with  the  first  term  In  (13.1?)  and  the 
latter  may  be  written: 


N 


A,cp^(y) 


1=1         M 


.a 


J=l 


tr  '  J' 


=  0, 


(15.19) 


vhere 


^ijW=-r^[^:W'^5w]- 


(15.20) 


"t-PJ 


It  Is  again  advantageous  to  \ise  vector  notation.   Due  to  the  linear  indepen- 
dence of  the  (f,,    (13.19)  is  equivalent  to  the  equation 


f(x)  +  W(x)a(x)  =  0:, 
from  which  it  follows  that 

a(x)  =  -  W"'^(x)f(x) 


and 


-1/ 


K(x,y)  =  (a(x),cp(y))  =  -  (W  ■'(x)f  (x),cp(y) ) 


(15.21) 

(15.22) 
(15.23) 


The  basic  property  of  the  elements  of  the  matrix  W(x)  is  similar  to  the  pro- 
perty (12.10)  of  the  elements  of  the  matrix  V(x): 


S  w,j(-)  -  ^(xXp^C.)  ■ 


(13.2't) 


This  relation  allows  one  to  again  obtain  an  elegant  expression  for  the  varia- 
tion in  the  potential  A  q(x)  and  the  solution  cpp(x,k): 


A  q(x)  «  -  2^  IQ  ||w(x)  II, 


(15.25) 


dx 
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cpgC^A) 


W(x)  f(x) 

3(x,k)     cp^(x,k) 

l|w(x)|| 


(15.26) 


vhere 


[q)P(x);cp,(x,k)]  PR 

p   (x,k)   =    '-J  ^  g      -^  =     /  cp?(t)cpT(tA)<it. 

'^  B .  +  k  J      J         -^ 

J  o 


(13.27) 


The  result   (13.25)   is   called  Bargmaim's   fonnula. 

In  Sec.    15^    It  will  be  shovn  that  cpp(x,k)    actually  is  a  solution  of 
(1.1)   -with- the  potential  q^(x)   =   q,  (x)   +  Aq.(x).      The  eisymptotlc   form  of 
cp   (x^k)    can  be  deduced  directly  from  the  explicit  form  of  (13.26)    Just   as 
vas   done  In  the  preceding  section.      It  Is  given  by 


cp2(x,k)  =  ^ 


21k 


where 


M2(-k)e^^  -  M2(k)e-^^ 


N 


+  0(1), 


Mgdc)   =  M^(k)        n       FTir  ^ 


(13.28) 


(13.29) 


and  this  proves  the  correctness  of  o\ir  solution. 

The  behavior  of   A  q.(x)    as  x  -►  0  or  x  -►  o»     can  easily  be   analyzed.      For 
X  -*  0  ,    A  q(x)   decreases  exponentially  as  before,   however,    as  x  -*  0      A  q(x) 
does  not  tend  to  zero.      Let  us   consider  one  factor  for  the  sake  of  simplicity. 
On  the  basis   of  (13.24) 


W(x)   =  W(0)  +     j      f^(x,ia)cpj_(x,i3)dx. 


(15.30) 
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2    2 
But  from  (l5.20),  It  follows  that  W(0)  =  M  (ia)/(a  -  (3  )  and  this  Implies 


M  (la)    /  2  \ 

W(x)  =  -| 2   1  +  (a  -  3  )  ^  [l  +  o(l)]    .  (13.31) 

a  -  P   \  / 


A(i(0)  =  -2(a^  -  p^).  (13.32) 


The  formulas  cited  can  be  used  to  determine  the  potential  approximately  for 
given  S(k) .   Any  S(k)  may  be  approximated  by  a  rational  function  S  (k) .   The 
corresponding  approximate  function  Mp(k)  will  be  given  by  a  product 


^   k  +  la 


M  (k)  =  n   i  .  (13.33) 

^  ii  ^^"-h 

The  potential  constructed  from  this  function  will  be  given  in  terms  of 
trigonometric  functions.   However^  In  general,  its  asymptotic  representation 
differs  from  that  of  the  exact  solution.   But  in  some  average  Interval,  the 
behavior  of  the  potential  will  be  described  by  the  approximate  solution  rather 
well. 

As  an  example,  take  the  function 

The  corresponding  phase  is  given  by 

eot  Ti(k)  =  ^2£_  l  +  _i.^k,  (13.35) 

'^'p-aka-p 

and  there  Is  no  discrete  spectrum.  For  the  potential  q(x)  one  obtains 

q(x)  =  2  P^(P^  -  ^  ) g  .  (13.36) 

(p  cosh  px  +  a  slnh  px) 

The  expression  for  cp(x,k)  is  more  Involved  and  will  not  be  given. 
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1^)-.   The  case  I  >  0. 


Equation  (l.l),  considered  up  to  this  time^  Is  a  particular  case  of 


(0 


y"  + 


^,(x)  .  iii-±JI) 


y  =  s  y^ 


(1^.1) 


which  arises  vhen  variahles   are  separated  in  the  three-dimensional  Schroe- 
dinger  equation  with  a  spherically-symmetric  potential  q(x)2 


Au  +  qu  =  8  u 


(llf.2) 


In  this  section,  all  the  properties  obtained  in  the  preceding  sections  for 

(I) 
the  operator  L  will  be  carried  over  to  the  operator  L   .   Let 


CP^^)(X,8)  =    1    J,(3X), 


l+l 


(1^.5) 


f^^^(x,s)  =  {±)^\[^\bx),  (111. if) 

where  y      ,    h^   (x)  are  spherical  Bessel  functions;  for  any  cylindrical  func- 
tion, the  corresponding  spherical  function  is  given  by: 


'^(x)  =  V^ 


2   Z  ^_^i(x)  {I   =  0,1.2,...). 


(1^.5) 


As  x-»0   or  x-*«>   these  functions  behave  in  the   following  way: 


(P^^^(X,8) 


l+l 


<P^^\x,s) 


X*»  00 


=    (21+1)11    [l+°(l)]^ 
=  -^  Sin  (bx  -  ~^+  o(l). 


J 


(1^.6) 
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(sx) 


(1^.7) 


f^*)(x,s) 


X*oo 


isx    ,^. 
=  e    +  0(1) 


They  play  roles  analogous  to  and  e    in  the  case  under  consideration. 

The  solutions  cp^   (x^s)  and  f^  "^(x^s)^  vhich  are  generalizations  of  cp(x,8) 
and  f(x^s),  are  determined  by  the  conditions 


11m  ^ j-^ — 9^  '(x,s)  =  1, 


x  -*  O    X 


l+l 


^.    -isx  Jt),        ,    . 
11m  e     f^  '(x,s)  =  1. 


(1^.8) 
(1^.9) 


X  -*  00 


The  integral  equations^    similar  to   (I.5)   aJid.  (I.6)    and  eq-ulvalent  to   (l^.l) 
with  the   conditions   (1^.8)    and  {lk.9),    are  given  byj 

X 

(p^^)(x,s)   =  cp[^'\x,s)   +   J    J^^)(s.x,t)q(t)cp^^)(t,s)dt,  (iJ+.lO) 

o 

00 

f(^)(x,s)  =  f^*^(x,s)   -    J    J^*)(s;x,t)q(t)f(^^(t,s)(it.  (1^^.11) 


Here, 


j(*)(s;x,t)  =  (is)* 


9^*^(x,s)f^*)(t,-s)   -9^*)(t,s)f(^*)(x,-s)].      (14.12) 


By  means  of  these  equations,  the  results  of  Sec.  1  carry  over  in  respect  to  the 
behavior  of  cp^   (x,8)  and  f^   (x,s)  In  the  complex  s-plane^  for  large  x,  etc.. 
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under  the  supposition  that  the  potential  q(x)  satisfies  the  condition 


/   x|  q(x)  |dx  <  00  - 


(1^-15) 


For  example,,  cp^   (x,s)  is  an  entire  function  of  a,    f^  {x,s)    is  analytic  in 
s  in  the  upper  halfplane  t  >  0,  and  for  these  solutions  the  ineq^ualities 


|cp^^^(x,s)|g  k( 


-2 V''\M^ 


1  +   S  X/ 


f^^^(x,s)|  ^  k(— ^— ^  )   e"''^   t  S  0, 


S  X 


( t) 
are  valid.     For  large    |s|,      f^    '^(x, s)   behaves   asymptotically  like 


(li^.ll+) 


(1^^.15) 


f^^^(x,s)  =  e^^^  +  o(l),   T  S  0. 


(Ii+.l6) 


For  real  s^  cp^   (x,k)  can  he  expressed  in  terms  of  f^   (x^k) 


-<"(-^)  =  afe(y' 


f^^^(x,k)M^^Vk)-(-l)V^^(x,-k)M^*^k) 


(1^.17) 


vhere 


M^'^^(s)  =  lim 


X  -  o  ±   (21-1) '.I 


(--)    f^^)(x.s) 


(11^.18) 


possesses  all  the  properties  of  the  function  M(s)  emmierated  in  lemma  1.6, 

vith  a  correction  in  the  formula  for  the  normalization  constant  C  .   In 

n 

the  present  situation^  C^  is  given  by 
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^i"  - 1 


CP^'^(x.iK 


'd:!c  = 


M(i/<    ) 

n 


n      ^ 


lim 
X  -►  o 


f^'^(x.l^ 


«,+! 


(li^.19) 


(I) 


,(0, 


As  before.    It  is   assumed  that  M^    '^(O)   /  0  since  the   condition  M^      (O)   =  0  is 
again  equivalent  to  instability  and  this   restriction  is   of  no   consequence. 
With  this  assumption^    the  solution  cp^      (x, O)   behaves  asymptotically  for  large 
X  like 


cp^^)(x,0)   =  Ax^^[l  +  o(l)]. 


(14.20) 


The  results  concerning  the  expansion  theorem,  the  existence  of  transforma- 
tion operators,  and  the  asymptotic  expansions  for  large  time  of  the  solution 
of  the  Schroedinger  equation  all  go  over  to  the  present  case.   Thus,  the 
completenesa  relation  for  the  eigenfunctions  is  given  by 

m 


2J  C^9^(x)cp^(y)  + 


n=l 


'  ''  ^    'dk  =  5(x-y). 


(li^.2l) 


Furthermore,  the  following  representations  for  cp^   (x,s)  and  f   (x^s)  hold: 


cp^^)(x,s)  =  ^[^\^,s)   +  J    K^^)(x,t)q,[^^^(t,s)dt, 


(l4o22) 


f^^)(x,s)  =  f[^\^,^)   +   f   A^^)(x,t)f[^^^(t,s)dt.  (li^.25) 
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/"  1. 


The  scattering  operator  associated  with  L^   is  constructed  with  the  help 
of  the  function 


(14.21;) 


This  same  function  occurs  in  the   asymptotic  representation  of  the  normalized 


,(0 


elgenfunction  t^    Hx,k)   =   (ik)  V      (xA)/m^      (k) : 


UO.      ,^/J0, 


^i'^^^(x,k) 


X-»  cc 


1 

2ik 


S(^)(k)e^^ 


,    ^.t     -ikx 
(-1)      e 


+  o(l).  (1^.25) 


As  "before^  these  restilts  allow  one  to  solve  the  inverse  problem^  i.e. 
the  reconstruction  of  L  from  its  S-function.  Thus  the  kernel  K  {x,y) 
satisfies  the  equation 

X 

K^^^(x,y)   +   n^^\x,Y)   +    J    K^^^(x;t)n^^^(t,y)dt  =0,      x  >  y,  (lif.26) 


where 


m 


n=l 


(^)^,.  v^v2(^l) 


I     f    q)^')(x,k)      ,^.       ^,. cp^'^(y,k)k 


dk. 


M^^(k)M^^(-k) 


(1^.27) 


By  a  repetition  of  the  earlier  argrmient,  this  equation  can  be  shown  to  have 

a  unique  solution.  Also,  the  completeness  relation  (ll+.2l)  for  Cp^  '(x,k), 

ft) 
constructed  from  the  solution  K^   (x^y)  'by  means  of  (lJ+.22),  can  be  proven, 

and  a  differential  equation  of  type  (l^t-.H)  can  be  derived  for  cp^   (x^k).  How- 

ever,  tracing  the  relationships  between  S^  '^(k)  and  the  potential  q(x)  using 
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this  equation  or  the  analog  of  the  Marchenko  equation  turns  out  to  be  diffi- 
cult.  This  is  due  to  the  fact  that  the  kernel  of  equations  of  type  (lij-.26) 
is  constructed  using  Bessel  functions  for  vhich  no  simple  addition  formula 
exists  as  for  the  trigonometric  functions.   Of  course^  conditions  like  (6.7) 
and  (6.9)  for  the  S-function  are  still  valid  in  the  present  case.   As  to  the 

analog  of  condition  (6.8);,  it  is  not  clear  beforehand  that  formulating  it  in 

( I) 
terms  of  the  Fourier  transform  of  S^   (k)  is  convenient.   One  might  think  that 

the  Fourier  transform  arose  in  a  natural  way  in  the  case  i  =  0  because  ve 

dealt  with  trigonometric  functions.  All  the  more  suiprising  is  the  fact  that 

(I) 
the  behavior  of  the  S-function  for  L^   turns  out  to  be  no  different  than 

that  of  the  S-function  for  L   .   We  have  in  particular 

Theorem  1^.1.   If  S(k)  is  the  S-function  for  the  operator  L    associ- 

ated  with  the  differential  equation  (l^.l)  with  a  potential  q(x),  then  it  is 

also  the  S-function  of  an  operator  L    with  arbitrary  m  =  0,1,2,..., t+1, . . . 

(m) 
where  the  corresponding  potential  q   (x)  behaves  like  the  potential  q(x)  as 

X  -►  0  and  x  -►  »  . 

The  proof  of  this  theorem  will  be  given  in  the  following  section.   Also, 

the  sense  in  which  the  behavior  of  the  potentials  q(x)  and  q   (x)  is  analogous 

will  be  made  more  precise  there. 


15.   Transformation  of  Sturm-Liouville  type  equations. 

In  the  preceding  sections,  it  was  repeatedly  mentioned  that  a  series  of 
assertions  would  be  proven  in  Sec.  I5  on  the  basis  of  a  certain  general  method. 
The  essence  of  this  method  will  first  be  presented  from  which  the  assertions 
will  then  follow. 
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Let  y  (x)  te  some  partlculaa:  solution  of 

-y"  +  q(x)y  =  Xy  (15.I) 

for  X  =  X     vhlch  does  not  vaMsh  In  the  neighborhood  of  the  point  x  =  a. 
Consider  the  expression 

[y(x,X);y  (x)] 

y  (X,X)  =  ^-—  ,  (15.2) 

(x-x,)y,(x) 

vhere  y(x,X)   Is   an  arbitrary  solution  of  (15«1)    and    [cpjt]    Is   the  Wronsklan 

[Cp;t]    =  CpKx).t(x)     -   Cp(x);|,«(x)     .  •  (15.5) 

Lemma  I5.I.      The   function  y^(x,X)   Is  a  solution  of  (l^.l)  vlth  a  potential 
q  (x)   =  q.(x)   +   Aq(x)   vhere 

''  o  dx 

To  prove  this,  we  note  that  any  two  arbitrary  solutions  of  (l5«l)  satisfy 

the  equality; 

[y(x,X  );y(x,X2)]    ,  ,    ,  ,  ,  , 

^    ,^_,^     =  y(x,.^)y(x,X2).  (15.5) 

Therefore, 

[y(xA))y„(x)] 

y'(x,X)  =  y(x,X)  -  — ~ y'(x)  =  y(x,X)  -  y  (x,X)v(x),         (15.6) 

'  (x-x^W 

where  the  function  v(x)  =  y'(x)/y  (x)  is  a  solution  of  the  Riccatl  equation 

v'(x)  +  v^(x)  =  q(x)  -  X^.  (15.7) 

Differentiating  (15 -6)  once  and  using  (15.?)^  we  find 
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ylix,X)   =  y'(x,X)  -  y|(x,X)v(x)  -  y^(x,X)v'(x)  = 

y;(x)       2 

=  y'(x,X)  -  y(x,X)  —j^  +   y  (x,X)v  (x)  -  y  (x,X)v'(x)  = 

=  - y  (x)     +  yi(x,X)[v2(x)  +  v'(x)J  -  2y^(x,X)v'(x) 

=  R  -  X  +  v^(x)  +  v'(x)  -  2v'(x)1y^(x,X). 


i.e. 

-  yJ[(x,X)  +  q^(x)y^(x,X)  =  Xy^(x,X).  (15.8) 

This  proves  the  lenma. 

Formula  (I5.2)  makes  sense  for  X  ^  X   .      If  X  =  X  .  one  of  the  solutions 
^  "^   '  '   o  o 

of  the-  transformed  ecjuatlon  Is  the  function 

As  a  second  linearly  independent  solution,  one  may  take 

X  X 

Of  course,  the  solutions  (15.9)  and  (15.IO)  can  be  deduced  from  (l5-2)  by  a 
limit  process.   ConTersely,  the  function  y(x,X)  is  expressible  in  terms  of 
solutions  of  the  transformed  equation.   F05:  this  it  suffices  to  note  that 
(15.9)  implies 

and,  therefore,  (15.6)  can  be  written  in  the  formi 
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y(.A)   =  y[UA)   -  y,(x,x)  ji^^  =  ^'  .^^(.r  -'  •     (15.12) 

This  expression  defines  the  transformation  inverse  to  (15.2). 

Till  nov^  ve  have  considered  the  transformation  (15.2)  in  the  neighbor- 
hood of  X  =  a  where  the  mapping  solution  y  does  not  vanish.   We  vlll  noV 
describe  the  behavior  of  the  solution  of  the  transformed  equation  when  this 
condition  does  not  hold.  Having  in  mind  its  application  to  (l.l)  and  (l^J-.l)^ 
we  will  assijme  that  the  singular  point  occurs  at  x  =  0  and  that  In  the  neigh- 
borhood of  X  =  0,  the  given  potential  has  a  singularity  of  the  type 

g(x)  =   '^'l^K  of-^Y      e>0.  (15.15) 

This  requirement  is  somevhat  stronger  than  the  condition  Imposed  on  the  poten- 
tial in  Sec.  1  and  Sec.  1^_,  namely 

oo 

/  x|q(x)  [dx  <  00  ^  (15.1^) 

o 

but  It  slB^llflea  considerably  all  of  the  calculations.   For  this  situation, 

(15.1)  will  be  said  to  have  an  l-slngularity  at  x  =  0.   Two  types  of  solutions 

exist  in  the  neighborhood  of  x  =  Oi  one  regular 

y{x,X)   =   C(X)x^+^[l  +  0(x^)],  (15.15) 

and  the  other  irregular 

z(x,x)  =^  [i  +  q(x^)].  (15.16) 

X 

Evidently,  all  reguleir  solutions  only  differ  by  a  factor. 

Consider  now  the  transformation  generated  MBlng,   some  regular  solution 
y(x,X  ) .  The  term 
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Aq(x)  =  -  2  ^  in  x^+^[l  +  0(x^)]  =  ^^-^^  +  ^(-^  ^^^-l?) 

dx  X        \x   / 

Is  added  to  the  potential  q(x)^  so  that  the  equation  has  an  -t-slngularlty  at 
x  =  0.  In  deriving  (15.I7),  the  asymptotic  representation  is  differentiated 
but  this  can  be  rigorously  justified. 

Let  us  observe  vhat  happens  to  a  regular  and  irregular  solution  under  our  . 
transformation.   By  virtue  of  (l5.o)  and  (I5.I5) 

[y(x,x)jy(x,x  )]    p 

~=    J     7{t,X)Y(t,\^)dt-  (15.18) 


X  -  X 
o 


There fore J 

X 

f    x2(^l)[l  +  0(x^)]dx 
y^(x,X)  =  C(X)  ^ ^ =  2!^  x^2[-^  ^  Q(^6)j^     ^^5^^5) 

X 

Furthermore, 

[z(x,X))y(x,X^)]    =   (2t  +  1)D(X)CCX^)[1  +  0(x^)]  (15-20) 

so  that 

z^(x,X)  =   ("^  ^^f  ^^^    [1  -H  0(x^)]  .  (15.21) 

X 

In  the  folloving,  (l5-2)  will  be  used  to  transform  a  regular  solution  and  (I5.I2) 
to  transform  an   irregular  one.   It  has  thus  been  shown  that  transforming  with  a 
regular  solution  increases  the  t-singularlty  at  the  origin  by  unity.  Moreover, 
regtilar  and  irregular  go  into  regular  and  irregular  solutions  respectively.   It 
is  easily  seen  that  transforming  with  an  irregular  solution  lowers  the  l-singu- 
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larity  of  the  equation  by  unity  but  again  preserving  the  regularity  property 
of  the  respective  solutions. 

ThuB^  to  obtain  an  equation  having  the  same  ^-singularity  as  the  given 
equation^  two  successive  transformations  must  be  performed^ the  first  using 
a  regular  solution  of  the  given  equation  and  the  second  then  using  an  irre- 
gular solution  of  the  transformed  equation.   It  turns  out  that  the  transform- 
ations constructed  in  Sec.  12  and  Sec.  15  hy  means  of  the  Gelfand-Levitan 
equation  could  be  deduced  in  this  way.  We  vlll  Illustrate  this  by  transform- 
ing tvo  operators  whose  spectra  differ  by  a  single  eigenvalue.   Consider  an 
equation  of  type  (l.l)  with  a  potential  satisfying  condition  (1.2)  and  assiune 

that  none  of  the  points  in  the  discrete  spectrum  of  the  associated  operator 

2 
L  lie  to  the  left  of  X  =  -  p  .   Take  some  p  >  p  .   Under  the  given  conditions^ 

(p  (x)  =  cp(x,ip)  vanishes  only  at  x  =  0.   Perform  the  transformation  (l5o2) 

using  this  solution.   Then  the  transformed  equation  will  have  an  ^-singularity 

at  X  =  0  with  ^  =  1.   The  function 


f^  (x)  =  — 7— r 


1  + 

o 


A. 

C  J   cPo(t)(it 


(15.22) 


will  be   an  irregular  solution  of  the  transformed  equation.      Performing  a  trans- 
formation with  this  solution,   we  arrive  at  an  equation  which  has  no  singularity 
at  x  =  0.      Let  us   calculate  the  potential  and  solution  cpp(x,k)    associated  with 
this  equation.     To  do  this,    it  is  necessary  to  combine  the  fonmilaa 

rcp(x,k);cp.  (x)|  ,  p 

^l(-^^)   =  T-2 2rTT  =  ^   J     ^(^'^)^o(*)^^'  ^^5.25) 

(k     +  p    )cPo(x)        ^o^         ^^ 

[cp,(xA);t,(x)] 
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,2  2 

Aq(x)   =   -   2~lJX  cp^(x)    -   2  ^  In  ^^^(x).  (15-25) 

dx  dx 

TMs  leads  to  the   follovlng  res\ilts! 

cp^(x)  p 

cp^Cx.k)   =  q>(xA) 2_ /     cp(t,k)9o(*)^^^  (15.26) 

1  +  C     I     cp^(t)dt     ° 
o 


Q(t)dt  ]     .  (15.2?) 


These  formulas  agree  vlth  (12.1^4-)  and  (12.15).  Ihus  by  a  direct  verification 
it  has  been  shown  that  the  function  cpp(x^k)  determined  by  these  formulas  is  a 
solution  of  (1.1)  with  the  potential  given  by  (l5.25).4-^ 

To  Justify  the  unproven  assertion  of  Sec.  13,    it  Is  necessary  to  perform 
the  following:  Transform  the  given  equation  using  the  solution  cp(x^l3).   Then 
transform  the  resulting  equation  using  the  solution  f  (x^ia)  obtained  from 
the  solution  f(x,ia)  of  the  given  equation  by  the  first  transformation.   By 
combining  the  formulas 

[cp(x,k);cp(x,ip)]      ^     r   ,    V  ,    X  /    .ox 

q,.(x,k)=— p P =  ^rTTBT  /  cp(t,k)q)(t,ip)dt,  (15.28) 

^        (k^  +  r)cp(x,13)   "P^^'^P^  i 


[f(x,la)jcp(x,lp)]     ,  X 
(p  -  a  )cp(x,13)    ^^  '    ^^ 

[cp^(x,k);f,(x.la)] 
^2^^'^)  =  f.^(x,ia) '  ^^^'^^ 
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Aq(x)  =  -  2  — ^  liicp(x,l3)  -  2  — ^  In  f  (x,ia). 
dx  dx 


(15.51) 


one  finds  that  the  function 


cp^Cx^k)  =  cp(x,k) 


f(x,ia)[cp(x,k)jcp(x,ip)] 


W(x) 


,  2   „2 


(15.52) 


is  the  solution  of  (l.l)  with  the  potential 


qgCx)  =  q(x)  -  2  ^  In  W(x)  =  q(x)  -  2  ~  In  -^ g- 


dT 
dx" 


dx 


2     [f(x,ia);cp(x,i3)] 

2 


3  -  a 


(15.55) 


This  proves  the  statement  of  Sec.  15. 

These  properties  of  a  transformation  vlll  now  be  used  to  prove  theorem 
1^.1.  A  transfonnation  using  any  regular  solution  changes  an  ^-singularity 
only  at  x  =  0.   The  behavior  of  the  potential  increment  Aq(x)  as  x  -♦  » 
may  be  different  depending  on  the  location  of  the  parameter  X  in  the  complex 
plane.  For  example,  in  the  cases  considered  till  nov,  the  potential  increment 
decreased  exponentially  as  x  -►  »  .  However,  there  exist  solutions  vhich 
change  the  singularity  of  the  equation  in  the  same  way  both  for  x  -»  0  and 
X  -*  00  .   Such  is  the  solution  of  equation  (15.1)  with  X  =  0.  If  it  Is  assumed 
that 

,        6  >0,  (15.5^) 


q(x)  =  ^ — 2 +  0' 

X 


\; 


2+5 


as  X -*  00  ,  then  using  equations  like(l4.10)  and  (l4.ll),  one  can  show  that 
the  regular  solution  has  the  asymptotic  representation 


y(x,0)  =  Cx 


l+l 


1  +  0 


(15.35) 
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as  X -♦  00  axid  among  the  irregular  solutions  there  exists  a  z(x^O)  such  that 


z(x,0)   = 


'*\7 


(15.56) 


as  X  -♦  00   .     Under  transformations  tising  these  solutions. 


Aq(x)   =   -   2 


dx  X  \x       / 


(15.57) 


21 


Aq(x)   =   -   2  ^  In   z(x,0)   =   -  -^  +  0 


dx 


2+6 


(15.58) 


Thus  they  have  the  requisite  property  of  changing  the  singularity  of  the 
equation  equally  for  x  -»•  0  and  x  -»  »  .  Let  us  see  hov  the  S-functions 
change  tinder  transformations  using  these  solutions.   It  is  easily  verified 
that  the  fiinctlons 


^     ^""'^^   -  -iky(x,0) 


-1)/   ^^     t^il!i±il!!i±^ 

-ikz(x,0) 


f(^-^)(x,k)  = 


(15.59) 
(15.^0) 


have  the  asymptotic  behavior 


Jt+1),     ,  V  ikx  /^v 

f^  (x^k)  =  e         +  0(1) 


(I5.i^l) 


as  X  -+  00  .      Thus  they  are  solutions   of  the  transformed  equations  of  type 
f'    '(x^k).      The  function  M(k)^    detennining  the  S-function,    arises  in  the 


asymptotic  formula  for  solutions   of  type 


f^^^(x,k) 


aB  X  -►  0 
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f^^\x,k) 


x-o  (kx) 


(I5.i^2) 


Hovever^   ty  virtue  of  ( 15.20)^   ve  find 


r^^^^(x,k) 


.    (21+1111  ,U1^(^) 
x-o       (kx)^"*"^ 


1  +  O(x^) 


(15.^3) 


f^^-^\x,k) 


x-»o 


.  (^^-^)i^  1^-Vk) 

(kx)^-l 


1  +  O(x^) 


(15.^^) 


Thus,  it  has  been  shown  that  under  transformations  using  the  solutions  y(x, O) 
and  z(x, O),  M  and,  consecLuently,  the  S-function  does  not  change.  This  asser- 
tion together  vlth  formulas  (I5.I7),  (15.57)  a^^d.  (15.58)  for  the  asymptotic 
behavior  of  the  potential  as  x  -*  0  and  x  -»  <»  ,  proves  theorem  l4.1. 
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APPENDIX 
(Comments  and  notes  on  the  literature) 

1..  The  proofs  of  lemmas  1.1  -  l.J  and  1.5  are  found  in  Levinson's 
paper  [6] .  Certain  statements  are  proven  in  papers  [27]  and  [28] .  The 
case  M(0)   =  0  is  detailedly  treated  by  Marchenko  and  Agranovich   [29],    [50] . 

2.      The   completeness  theorem  for  the  eigenfunctions   of  the  operator  L 
in  the  form  (2.6)   is  proven  by  LevlnBon    [6]    in  the   absence  of  a  discrete 
spectrum.      The  general   case  is   considered  by  Jost  and  Kohn    [9] • 

5.     Many  papers  have  been  devoted  to  the  question  of  hov  the  solution 
of  the  time  dependent  Schroedinger  equation  behaves  for  large    |t|   both  in 

physics   and  mathematics.     A  nonrigorous  proof  of  the  existence  of  limits 

iLt  "^^o"*^ 
of  the  operator  U(0,  t)3se       e  as  t  -^  ±  «>  ,    typical  of  physics  papers, 

is   cited  for  example,    in  the  survey  of  GeHmaxin  and  Cxoldberger   [51]  •     From 
mathematical  vork,    it  is  necessaiy  to  mention  first  of  all  the   articles  of 
Friedrlchs    [25],   who  proved  the  existence  of  the  limits  U(0,±'»).      He   also 
shoved  that  the  limiting  operators   are  unitary  for  a  vide   class   of  unper- 
turbed operators  L     on  the  assumption  that   the  perturbation  operator  V  is 
small,     A  formal  presentation  of  his  method  is   cited  in  the  paper  of  Moses 
[52].      Cook    [55J   proved  the  existence  of  the  limits  u(0,±«>)    for  the  three- 
dimensional  operator  -  Au  +  q(x)u  in  all  space  assumi ng  only  that,  q(x)   is 
square  Integrable  in  all  space.      Hovever,    he  did  not  study  the   question  of 
vhether  the  operator  S  =  U(0,«')-'<-U(0, -«>)   is  iinitary.      The  restriction  that 
the  perturbation  operator  V  be  small  is   removed  in  the  paper  of  Ladlzhenskaya 
and  Faddeyev   [54]   using  the  formalism  of  Friedrlchs. 
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Theorem  3.I  does  not  follov  from  the  results  of  these  papers  for  the 
restrictions  we  have  imposed  on  the  potential  q(x).   The  elementary  proof 
cited  makes  use  of  the  concrete  properties  of  the  example  under  considera- 
tion and  does  not  carry  over  to  other  problems. 

^.   Povzner  [55]  aJ^d.  Levitan  [56]  first  obtained  and  used  the  repre- 
sentation (^.5)  for  cp(x,k).   Formula  (h,2)   and  f(x,k)  was  deduced  by  Levin 
[37]  Slid   our  first  derivation  repeats  his  argument.   The  method  of  deriving- 
the  integral  equations  (^.4)  and  (4.5)  as  well  as  the  inequalities  (4.7) 
and  (4.8)  is  due  to  Agranovich  and  Marchenko  [29]^  [50] . 

The  theorem  of  Titchmarsh  mentioned  can  be  formulated  in  the  follow- 
ing way:  A  necessary  and  sufficient  condition  for  ^(x)  to  be  the  limit  of 
some  function  J(z)  =  J(x  +  iy)^  which  Is  analytic  in  the  upper  halfplane  and 
such  that 

op 


f     1  $(x  +  iy)|2dx=  0{e-^^), 


la  that 


(p(t)  =  ^  5(x)e"^^dx  =0^    t  <  k. 

—  so 

5-   The  general  concept  of  transformation  operator^  as  already  noted^ 
was  developed  by  Friedrichs  [22],  [23]'   Some  of  the  notation  and  the  proof 
of  theorem  5.2  was  taken  from  the  articles  of  Kay  ajid  Moses  [l8]^  [2o]  who 
applied  Friedrichs'  method  in  solving  inverse  problems. 

From  the  formula  (5.I5)  for  the  S-functlon,  it  follows  that  S(k)^  in 
general,  cannot  be  continued  into  the  complex  k-plane.   Thus,  Heisenberg's 
supposition  that  the  discrete  energy  levels  in  the  example  in  ques^ion  might 
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be  determined  by  the  analytic   continiiation  of  the  S-matrlx  is  not  Justified, 
This  fact  was  noted  by  Jost    [27] . 

The  differential  equation  (5.25)  and  the  condition  (5.2^)  is  the  lead- 
off  point  for  the  proof  of  the  existence  of  the  kernel  K(x^y)  in  the  paper 
of  Gelfajid  and  Levitan  [12]  .  One  easily  obtains  th&  integral  equation  (^.^) 
from  this  equation.  Chudov  [6^]  proposed  using  the  nonlinear  equation  got- 
ten from  (5.25)  by  replacing  q(x)  by  2dK(x^x)/dx  to  solve  the  Inverse  prob- 
lem. Giving  the  S-function  for  large  x  provides  Cauchy  data  for  thi©  equa- 
tion. 

6_.      The  method  of  relating  W(k)   and  S(k)    on  the  basis  of  the  Wiener- 
Levi  theorem  vas  due  to  Krein    [l6]j    [58] •      The  Wiener-Levl  theorem  can  be 
stated  In  the   following  way.      Let  the  function  J(z)  be  analytic  in  a  region 
D  and  then  let  F(X)   be  so   chosen  that  the   curve  z  =  F(X)    (-  «  g  X  ^  <»)   lies 
inside  D.      If  F(X)   Is   representable  in  the  form 


f(x)  =  C  +     /     f(t) 


e       dt 


where  f(t)   is   absolutely  Integrable^   then  ^(F(X))    also  possesses  this  pro- 
perty. 

Formula  {6.9)  was   deduced  by  Levlnson    [6]    and  is   called  Levlnson's 
foiraula. 

7.  The   cited  relationship  between  the  kernels  K(x,y)   and  A(x,y)    does 
not  appear  in  the  literature. 

8.  The  first  derivation  of  (8.5)   is  taken  from  the  paper  of  Kay  ajid 
Moses    [18].     The  derivation  of  the  general  equation  (8.l4)   repeats  the  argu- 
ments of  Gelf and  and  Levitan   [l2] . 
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9.  The  existence  proof  for  (8.1^)  is  taken  from  the  paper  of  Jost 
and  Kohn  [ij]  and  to  a  great  extent  follows  the  reasoning  of  Gelfand  and 
Levi tan.  The  subsequent  presentation  with  certain  modifications  repro- 
duces the  arguments  of  Kay  and  Moses  [l8] . 

10.  The  analysis  of  the  properties  of  the  potential  q(x)  is  taken 
from  the  monograph  of  Agranovlch  and  Marchenko  [jo] .   The  various  rela- 
tionships between  W(k)  or  S(k)  and  q.(x)  were  obtained  by  Neiihaus  [39]  • 
Friedman  [4o] ;,  Jost  [^l] ,   and  Newton  [^2]  .  However  the  mor^  general  re- 
sults follow  from  (lO.ll)  and  (l0.2^). 

11.  Krein's  methods  are  given  in  a  series  of  articles  [^5]^  [^^] ^ 
[l6j^  [45]  in  DAN  SSSR  (see  also  the  lectures  presented  at  MGU  in  I956- 
1957) ♦   Only  certain  of  his  results  are  cited  in  the  survey.   The  system 
of  differential  equations  (II.I9)  is  the  starting  point  of  Krein's  methods, 

12.  Formula  (l2.l4)  for  the  Increment  in  the  potential  was  obtained 
by  Jost  and  Kohn  [15]  .   The  simplest  formula  for  the  solution  cp(x,k)  of 
the  type  (12. I5)  was  due  to  Kreln  [l6j  (for  the  case  m  =  1). 

The  portion  of  theorem  12.1  relating  to  necessary  and  sufficient 
conditions  is  due  to  Marchenko  and  Agranovlch  [29] ,    [50]  . 

15.  In  solving  (15.5)^  we  have  followed  the  paper  of  Fulton  and 
Newton  [i^•6]  who  refer  to  the  work  of  Bargmann  as  the  source  of  the  method 
uBed.  Expression  (15.25)  is  called  Bargmann 's  formula.  Formulas  for  cp(x,k) 
such  as  (15.26)  are  cited  by  Thelss  [^T] . 

Another  approach  to  the  problem  is  developed  by  Kreln  for  R,(k)  =  1. 
The  results  obtained  are  decisively  formulated  in  [48j ,  the  formulas  being 
simpler  than  (15.25)  and  (15.26).  However,  their  generalizations  to  the 
case  M^  (k)  ^  1  are  not  deduced. 
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1^.      The  basic  properties   of  solutions  of  (l^.l)   for  I  >0  are  obtained 
by  Levinson   [6] ,   Jost   and  Kohn    [ij] ,    and  Nevton    [^9] .      The  papers  of  Stashev- 
skaya    [50]    and  Voli    [51]   are  devoted  to   carrying  over  the  results   of  Gelfand 
and  Levltan  to  equations  vith  a  singularity  at  x  =   0.      Theorem  1^.1  is  due 
to  Marchenko   (it  was  presented  at  the  Aprils    1956  meeting  of  the  Kharkov 
mathematical  society). 

15.     A  transformation  of  the  type   (15.2)  -was  applied  for  the  first  time 
by  Crum   [52]^   vho  used  it  to  change  a  differential  operator  over  a  finite 
interval  into  an  operator  having  one  less  eigenvalue  than  the  original  opera- 
tor.    Krein  extended  Crum's  method  and  applied  the   results  to  get  a  complete 
characterization  of  the  spectral  function  of  an  equation  vith  the  singular- 
Ity  l{l  +  l)/x     in  the  potential  at  x  =  0.     Marchenko  made  use  of  an  analo- 
gous transformation  to  analyze  the  relationship  between  the  S-function  and 
a  potential  given  by 

q(x)   A  q   (k)    -f  iiiiil  ,      r   x^*|q^(x)ldx<oo    . 

o 

The  presentation  in  the  survey  differs  somewhat  from  the  methods  of  the  above- 
mentioned  authors. 

It  Is  interesting  that  inasmuch  as  we  verify  formulas   {l2..lk) ,    (12. I5) 
and  (15.25),    (13.26)  by  algebi*aical  means  without  reference  to  the  theory  of 
the  general  Gelfand-Levitan  equation,    they  still  hold  for  complex  values  of 

the  paxameters   k  ,   C  ,  a.   and  p   .      The  associated  potential  is,    generally 

0 
speaking,   a  comgplex  function  with  a  singularity  of  the  type  m(m  +  l)/(x  -  x  ) 

at  those  points  where     ||v(x)  1|    =  0  or    |lw(x)  ||    =  0.     Here,  m  is  the  multi- 
plicity of  any  such  existing  zero.     This  fact  was  ittoted  by  Krein   pufi]    and 
Theiss    [if?]. 
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We  nov  briefly  consider  some  of  the  generalizations  of  the  problem 
investigated  in  our  survey.   By  analyzing  many  of  the  formulas  in  the  text^ 
one  sees  that  with  appropriate  changes  they  remain  valid  for  systems  of 
equations^  i.e.  for  the  matrix  equation  generalizing  (l.l): 


-  Y"  +  Q(x)Y  =  k^  . 

Here^  Q(x)  is  a  real  symmetric  matidx.   The  solutions  cp(x^k)^  f(x^k)  and 
the  functions  w(k)_,  M(k)^  and  S(k)  become  matrices  in  this  case.   Therefore 
it  is  necessary  to  pay  attention  to  the  order  of  factors  in  generalizing 
formulas  to  the  matrix  case.   The  matrix  M(a)  is  analytic  in  the  upper  half- 
plane  t  >  0  and  singular  at  those  points  corresponding  to  the  discrete  spec- 
trum.  The  matrices  W(k)  and  S(k)  are  related  to  it  by  the  formulas 

W(k)  =  M(k)"'S^'^(-k  )"■'",   S(k)  =  M(-k)M'^(k)"'^, 

M  (k)  being  the  transposed  matrix.  Similar  systems  were  studied  by  Jost 
and  Newton  [55]^  Krein  [^5]^  and  Agranovich  and  Marchehko  [29]^  [5oJ  •  A 
fundamental  difficulty  arises  in  carrying  over  the  discussion  of  Sec.  6  to 
the  matrix  case.   In  consequence  of  the  noncommutativity  of  the  matrices^ 
the  forcmilas  cited  there  no  longer  hold.  Finding  how  the  matrices  W(k) 
and  S(k)  are  related  reverts  to  studying  integral  equations  of  the  typet 


K(t)  =  F(t)  +  /  F(t  +  s)K(s)ds. 


o 

Marchenko  and  Agranovich  deduced  necessary  and  sufficient  conditions  on 
the  S-matrix  so  that  it  corresponds  to  a  matrix  potential  Q(x)  from  a 
given  class  making  use  of  analogous  integral  equations.  The  formulation 
of  conditions  directly  in  terms  of  the  S-matrix  still  remains  an  unsolved 
problem. 
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Nevton  [^9]  and  Agranovich  and  Marchenko  [5^]^  [5o]  considered  a 

o 
system  In  vhich  the  potential  has  the  singularity  I  (t  +  1)6  /x  j 

Agranovich  and  Marchenko  in  studying  such  a  system  reduced  it  to  a  regu- 
lar one  by  a  transformation  generalizing  the  kind  described  in  Sec.  15 • 
The  inverse  problem  for  a  system  has  mainly  been  treated  with  the 
aim  of  seeing  what  means  are  necessary  for  solving  the  inverse  problem 
for  the  Schroedlnger  equation 

2 
-  Au  +  q.(x)u  =  k  u 

in  nl  1    of  space  when  the  potential  decreases   in  all  directions.     However, 
this  problem  essentially  differs   from  those  treated  till  now.      In  fact, 
the  S-matrix  in  this   case  is  determined  by  the  so-called  scattering  ampli- 
tude f(k;a, p)   depending  on  the  wave  number  k(0  S  k  <  «)   and  two  unit  vec- 
tors a  and  p.      Thus   the  S-matrix  depends   on  a  larger  number  of  parameters 
than  the  potential  q(x)   which  may  be  regarded  as  a  function  of  the  distance 
r     (0  S  r  <  00)    and  one  iini  t  vector.      In  this  sense,   the  problem  is  over- 
determined  and  It  is  necessary  to  look  for  nontrlvlal  properties  of  the 
S-matrix  which  would  decrease  the  amount  of  parameters   on  which  it  depends. 
The  simplest  problem  where  an  analogous   situation  occurs  is  in  recon- 
structing a  decreasing  potential  In  the  one -dimensional  Schroedlnger 

equation : 

2 
-y"   +  q(x)y  =ky       (-«><x<oo) 

from  Its  S-matrlx.     In  this  case,   the  S-matrix  is  a  2-by-2  matrix: 
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S        (k)     B        (k) 

\  /  a(k)      b(k) 

s(k)  =■ 

,b(k)      c(k) 


and  due  to  the   condition  of  unltariness,    is  determined  by  giving  three   real 
functions   of  k     (O  S  k  <  »).      The  potential  can  be   regarded  as  being  given 
by  tvo  real  functions   of  x     (o  ^  x  <  «). 

The  inverse  problem  for  this   case  vas   considered  by  Kay  and  Moses    [20]^ 
[21]    (as   an  example  illustrating  their  general  approach  to  the  inverse  prob- 
lem)   and  by  the   author    [55] .      In    [55]    it  is  shown,   that   an  additional  condi- 
tion on  the   S-matrix  follows  from  the   analyticlty  of  the   coefficient  b(k) 
in  the  upper  hali^lane  t  >  0.      This   condition  implies  that  every  S-matrix 
(and  potential)   is  determined  by  one  of  the   coefficients   a(k)   or  c(k)  which 
may  be   chosen  as   an  arbitrary  function.      The  reconstruction  of  the  equa- 
tion with  an  arbitrary  potential  over  the  entire  aria   from  its  spectral 
matrix  function  was  treated  by  Bloch   [65] • 

The  elements  of  the  S-matrix  are  also  analytic  in  the  three-dimensional 
case.      The  proof  of  this  fact  is  given  in  the  papers  of  Khuri    [56]    and  the 
author   [57]    in  connection  with  the  so-called  dispersion  relations. 

It  is   interesting  to  note  in  this   connection,    the  statement  of  the  three- 
dimensional  inverse  problem  as  proposed  by  Moses    [58]  :      The  potential  q_(x) 
has  to  be  determined  from  its  back  scattering  amplitude  g(k,a)   =  f(k)a,-a) 
where  a  is  a  vector  running  over  a  hemisphere.     These  two  pieces  of  data, 
namely  two  real  functions  of  k     (O  ^  k  <«)    and  a,   which  runs   over  a  hemi- 
sphere,   contain  as  many  parameters  as  does  the  potential.      It  is  very  probable 
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that  the  procedure  of  Moses  converges  for  sufficiently  small  g(k,a)  which 
in  other  respects  may  be  a  quite  arbitrary  function. 

Ihere  exists  a  series  of  papers  in  which  the  inverse  problem  is  solved 
for  relativistic  equations  when  the  latter  reduce  to  ordinary  differential 
equations.   The  equation  obtained  by  separating  variables  in  the  Klein- 
Gordon  equation  has  been  studied  by  Corinaldesi  [59] .   The  one-dimensional 
Dirac  equation  was  considered  by  Kay  and  Moses  [60J ,    Toll  and  Prats  [61] 
and  Verdi  [65] •   In  all  of  these  papers^  a  relationship  is  established  be- 
tween the  potential  and  asymptotic  phase,  the  latter  being  given  for  both 
positive  and  negative  energies.   This  data,  just  as  in  the  problems  described 
above,  depends  on  a  larger  number  of  parameters  than  does  the  potential.  A 
correct  formulation  of  the  problem  for  the  radial  relativistic  equation  till 
now  has  not  been  given. 
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